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THE VORTEX CHAMBER AUTOMATIC FLOW-CONTROL DEVICE 


SYNOPSIS 


Coefficients discharge for flow through horizontal orifices are greatly 
reduced when vortex present. Several flow control devices have been de- 
signed utilizing this phenomenon. This paper presents rational explanation 
the resulting coefficient variation from dynamic similarity. The results 
laboratory experiments are shown, giving means predicting orifice flow 
characteristics from the inlet velocity condition. 


INTRODUCTION 


The free vortex has long been significant phenomenon 
the design hydraulic structures. For the most part, however, has been 
regarded undesirable possibility, and its consideration has been limited 
the design expensive appurtenances intended eliminate, least re- 
duce, the circular movements the fluid. 
recent years, the properties the free vortex have been utilized suc- 
cessfully solve several flow control problems. made studies 
“counterflow brake” for use pumping installations. This device consists 
spiral shaped chamber with tangential and axial connections. The normal 
direction flow the center (axial) connection and out the peripheral 
(tangential) connection. the flow reversed due power failure, the re- 
sulting tangential inlet conditions cause vortex motion within the chamber 
and therefore greater resistance flow. This device thought type 
pipe fitting which may replace the check valve special cases. 


Note: Discussion open until May 1960. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. Paper 2310 
part the copyrighted Journal the Hydraulics Division, Proceedings the 
American Society Civil Engineers, Vol. 85, No. 12, December, 1959. 
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Stevens(5) made practical use the vortex phenomenon the design 
diversion structure for combined sewer system. small dam the sewer 
the downstream side the diversion structure diverted all the low sani- 
tary flow into orifice chamber. The radial approach the horizontal ori- 
fice resulted normal orifice head-discharge relations. The high flow rates 
the main sewer for storm condition would pass over the dam, causing the 
entrance the diversion chamber tangential and vigorous vortex 
form. The resultant reduced orifice coefficient discharge allowed only 
small portions the storm water pass into the interceptor. 

Both these control devices exhibit the advantage relatively simple con- 
struction and maintenance because they have moving parts. should 
expected that these ideas would find more popular use hydraulic design ex- 
cept for the fact that date very little information design value availa- 
ble concerning the characteristics orifice flow with significant vortex 
present. 


Theory 


The following signs and symbols are used this paper and have been listed 
here for easy reference. 


area orifice 

diameter boundary ring 

orifice coefficient discharge 

diameter orifice 

centrifugal force per unit volume 

inertia force per unit volume 

acceleration due gravity 

static head center orifice 

pressure 

Reynolds number 

tangential component velocity vector 
discharge component velocity vector 
velocity vector 

vortex number 


co-ordinates point the flow field 


angle the divergent jet makes with the horizontal 
circulation 

mass density 

dynamic viscosity 


has been shown dimensional analysis(5) that the variables affecting 
flow through horizontal orifice are known 
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then the discharge equation will take the form 


(1) 


the circulation for the free vortex flow and defined 


The drawdown the water surface resulting from this vortex motion (see 


Fig. 
(4) 


Figure 


Water Surface 


Orifice Plate 
Divergent Jet 
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The resulting water surface profile hyperbolic. Superimposing radial 
flow this vortex motion does not change the shape the curve except the 
regions immediately near the orifice, because both components the resul- 
tant velocity vary inversely with the radius. 

the earlier studies made the University Wisconsin, (5) was noted 
that for vortex conditions strong enough alter the orifice coefficient dis- 
charge, the angle which the velocity vector makes with the radius, such that 
the radial component negligible value computing the drawdown. 

studying the water surface profile, and computing the circulation, 
from Eq. (5), the coefficient discharge Eq. (2), was found 


between limits 0.8 and The upper limit amounts extrapo- 
lation, since data were accumulated beyond 2.7, approximately. 
might expected that deviations from this relationship would found 
higher values felt that greater deviations the actual profile from 
the theoretical should encountered this region because greater shear 
effects. 

laboratory data have been accumulated previously which would allow 
computation the vortex number, from the actual chamber flow pattern. 
Measurements the velocity profile various points within the orifice 
chamber using immersion type meters are made difficult because their dis- 
turbance the vortex motion. This fact leads legitimate criticism the 
design value the curves and equations which are based assumed 
velocity-surface profile relationship. 


(5) 


Dynamic Similarity 


studying the dynamic similarity flow through metering devices when 
the forces governing the particle motion are the surface forces represented 
and the viscous forces, which combine give accelerative affects indi- 
cated the inertia force, fj, the parameter which used determine simi- 
larity different flow cases the Reynold’s number, For flow through 
particular device such the pipeline orifice, the coefficient does not vary 
greatly with except the regions low Reynold’s numbers. 

This parameter described being the ratio the inertial forces the 
viscous forces. However, the division one vector another has known 
mathematical physical significance; hence, the relationship expressed 
the absolute scalar ratio the two forces. Besides insuring similarity 
force magnitudes, necessary insure similarity kinematic effects 
(velocity, acceleration) represented the flow pattern, which amounts 
accounting for the vector aspects the problem. Ordinarily this accom- 
plished specifying the geometric similarity with particular values the 
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ratio orifice diameter pipeline diameter. This permissible since for 
given boundary configuration, the flow pattern uniquely defined (if suf- 
ficient approach length, measured usually terms orifice diameter, 
specified order insure uniform flow the approach the orifice. This 
another geometric effect.) 

the case orifice flow with the possibility tangential components 
the velocity vectors, similarity cannot defined terms the scalar 


Reynold’s number and the boundary geometry indicated because for 


particular total acceleration indicated the Bernoulli equation solution, 


(8) 


there infinite number possibilities for the direction the vectors. 
This problem can handled breaking the total inertial force into com- 
ponents and specifying 


the ratio 


(11) 


will determine kinematic similarity, and along with B/D (geometric), and 
dynamic similarity will assured. This ratio will seen the square 
the vortex number Eq. (2) or, since the velocity which characteristic 
the inertial effects the discharge velocity, and the value which 
characteristic the centrifugal component (Eq. (3)), also the 
square Eq. (7). 


Experimental Methods and Results (Fig. 


The experimental work was done the Marquette University Hydraulics 
Laboratory and consisted tests made sheet metal tank feet square and 
feet high. Orifices were cut 2-foot diameter, gage brass plates and 
mounted the center the bottom the tank. Two inner boundary rings 
were used. Each was constructed gage sheet metal and had two diame- 
trically opposite openings inches wide which were attached parallel sheet 
metal plates forming the entrance channel, which were fixed position mak- 
ing arbitrarily chosen angle degrees (or degrees for some the 
tests) with the tangent the boundary ring. The larger ring was 3.5 feet 
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diameter and 2.5 feet high and the smaller ring was 2.0 feet diameter and 
2.5 feet high. Water was admitted the tank through pipes inches di- 
ameter each which branched near the floor the tank into perforated 
inch pipes which distributed the inflow uniformly around the outer boundary. 
Vortex profile measurements were taken check the earlier University 
Wisconsin with moving point gage mounted aluminum beam 
feet long with 2-1/2 inch 5/8 inch cross section. Discharge was 
measured 3-inch orifice 4-inch pipeline upstream from the test area. 
For the low-flow tests 2-inch nutating disc type water meter was used. 

the current series, test runs were made which the head was varied 
from 0.3 1.8 feet, varied from 0.85 2.6, and varied from 0.2 0.48. 

For the initial tests the larger boundary ring was used and the entrance 
channel vanes set degree angle with the boundary tangent. Tracers 
and small current meter were used check for uniform velocity within 
the entrance channels. significant variation velocity from the inside 
outside the channels was found; therefore, the average velocity was com- 
puted from the discharge and cross-sectional area, and this value was used 
further computations. Within the boundary ring the velocity was found 
essentially tangential (negligible radial components except immediately near 
the orifice) and its magnitude the same the entrance velocity. The change 
direction the fluid velocity vector had thus been accomplished with 
appreciable change momentum indicating negligible shear forces. For 
some the tests, both entrance channels were used, and for others one was 
sealed shut. The results either case were identical even though some ec- 
centricity the air core was noted for the one vane condition. 

The results for the degree approach angle exactly corresponded the 
degree angle. The velocity within the boundary (at 1.75 foot radius from 
the center the orifice) was found the same that within the entrance 
vanes, and this value was used computing the circulation from Eq. (3). 


From Entrance Velocity 


Profile Tests 


Orifice Coefficient, 
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The curve (Fig. shows the variation orifice coefficient with vortex 
number for these data computed from the entrance velocity. The data 
closely follow the straight line indicated the water surface profile measure- 
ments for vortex numbers less than 2.0. Above this value separate curve 
evidently established. This deviation probably occurs because the profile 
computations were based measurement made only one diameter away 
from the orifice, more closely approximating the true energy condition that 
region. High vortex numbers represent conditions high velocity near the 
orifice and therefore greater shear losses already noted. 

further check the normal acceleration relations, measurements 
the angle (Fig. the divergent jet were made. From Eq. (3) 


(because upon release, the circular motion becomes tangential) Substitution 
Eq. (7) gives 


(13) 
(14) 
Results these measurements are shown below. 
0.952 70° 
1.25 66.6° 64° 
1.70 57.3° 60° 
1.86 53.7° 52.5° 
2.24 44.5° 46° 
21.25 49° 


Figure 


Experimental Apparatus 
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Because the small clearance between the tank bottom and the floor the 
Hydraulic Laboratory, considerable difficulty was experienced making these 
measurements. complete series such tests has thus not yet been attempt- 
ed. Because the jet was not well defined (considerable wavering) the check 
indicated these data better than had been expected. 


CONC LUSIONS 


Prediction orifice coefficients for vortex flow can made the basis 
the velocity inlet the orifice chamber. The accuracy such prediction 
for chambers which are concentric with the orifice should depend the abili- 
estimate the inlet velocity. The curves shown should extended in- 
clude larger vortex numbers. this may necessary build ap- 
paratus which the head and inlet velocity can independently controlled 
(such model the Stevens combined sewer system which affords the ad- 

vantage for laboratory purposes by-passing part the total flow. 
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THE SETTLING PROPERTIES SUSPENSIONS 


ABSTRACT 


This paper contains the results analytical and experimental investi- 
gation the settling properties suspensions particles fluid. The use 
these properties predicting the sedimentation the particles outlined. 


INTRODUCTION 


many engineering problems necessary deal with flowing fluid 
which particles are suspended. The engineer commonly encounters such 
problems river channels, water and sewage clarification, reservoirs 
and the delta regions rivers. One the tasks associated with these 
problems determine how much material will settle out suspension and 
where that material will settle. 

The factors that affect the settling suspended particles can divided 
into two groups. Those the first group can called the conditions flow. 
They are the temperature and pressure the suspension, the velocity distri- 
bution the flow, and the nature the turbulence flow. The factors the 
second group can called the settling properties the suspension. 

The term “settling properties suspension” refers how the particles 
behave given set conditions flow. Thus, for some specified tempera- 
ture, pressure, fluid velocity and turbulence, the particles suspension 
will settle, flocculate diffused some manner. The manner will vary 
from suspension suspension, and must determined experimentally for 
each. 

Hence, the problem determining how material will settle out flowing 
suspension can divided into two parts. First the flow conditions must 
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measured, and second, the settling properties the suspension must mea- 
sured. Neither these two easy task, and considerable research and 
development necessary before either will done satisfactorily. 

This paper deals only with the settling properties suspension and their 
use calculating the removal the particles. Section contains general 
discussion the problem calculating the removal. This discussion points 
out how the settling properties are used and indicates the type experi- 
mental information needed. The experimental determination settling pro- 
perties discussed Section Some results already obtained are present- 
ed, and the experiments yet developed are outlined. Section this 
research related more practical engineering problems. 


Calculating the Removal Suspended Particles 


The Problem 


was stated above that one the tasks the engineer determine 
advance how the particles will settle out given flowing suspension. 
specific, consider suspension particles water which flowing 
open channel. The fluid velocity and turbulence level are low enough for the 
particles settle the bottom the channel. The problem calculate 
the amount settled material function distance along the channel. 

Such calculation involves the solution some form continuity equation 


expressing the conservation suspended matter any point the suspension. 


order derive the equation first necessary consider the factors af- 


fecting the concentration particles point. The first these factors 
the settling the particles. 


Settling the Suspended Particles 


The term velocity suspended particle” refers the velocity 
that the particle would have were settling perfectly still fluid. sus- 
pension may contain particles many settling velocities and the distribution 
these settling velocities important settling property the suspension. 
calculations related settling just reasonable characterize par- 
ticles these settling velocities their size, shape, and density. For, 
the settling velocity that indicates how soon the particle going re- 
moved from suspension settling the bottom. 

Let represent settling velocity and let positive downward. The ve- 
locities the particles suspension will range from zero the velocity 
the fastest particle. For numerical calculations and theoretical discussions 
convenient divide the range values into number increments, 
Aw. Then can refer class settling velocities between the values 
and Aw, and those velocities between and 2Aw. general, will 
indicate velocities between the values and iAw. The particles with 
velocity will called i-particles. 

The concentration i-particles any point the suspension will 
called fj, the units which are mass per unit volume suspension. This 
concentration can vary from point point and from time time. Letting 
and represent coordinates which locate point, and represent time, the 
concentration will function, The term indicates that fis 
function which varies varies. Hence, the frequency function for the 
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SUSPENSIONS 
settling velocities point x,y,z. this paper will called the settling ve- 
locity distribution x,y,z. 

any point the suspension, the flux i-particles due settling will 
given 


This the rate which i-particles pass through unit horizontal area 
point x,y,z, and time The rate which all particles pass through this area 


co 


definition, this total flux simply the product the local instantaneous 
mean settling velocity, w(x,y,z,t), and the local instantaneous particle concen- 
tration, That is, 


where 


Hence, the product describes the motion particles due settling. 


Movement Particle Due Fluid Motion 


any point the suspension, the particles also have motion due fluid 
flow the point. the flow without turbulence, this motion simply the 
resultant fluid velocity. Let this resultant represented the vector 
U(x,y,z,t). will have x,y, and z-components U(x,y,z,t), V(x,y,z,t) and 
respectively. 

When the flow suspension turbulent, the motion the particles 
usually considered two parts. First, the particles are considered being 
carried along the temporal mean fluid velocity, also called U(x,y,z,t). 
Superimposed this motion the diffusion particles the turbulent 
fluctuations fluid velocity. The decision what the main flow and 
what turbulence will depend upon the flow pattern for each individual case. 


order deal with the particle motion caused turbulence, common 


define coefficient diffusion. Let the coefficient for i-particles 
and direction. Then the flux i-particles through unit area normal the 
direction can given 


Since there may also ejy and eyz, the coefficient will called when 
specific direction implied. 

The coefficient should not confused with the coefficient for diffusion 
transfer momentum between neighboring elements fluid. For the sake 
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differentiation the latter will called reasonable assumption, often 
made, that equal For detailed discussion this point, the 
reader referred the experimental and analytical work and 
Their work indicates that not only differs but that the 
manner which differs varies with For coarse sand large) flowing 
water, was smaller than e,, while for fine sand small) was larger. 
Nevertheless, the order magnitude the difference was such that may 
reasonable assume that engineering calculations. With this 
assumption there need for the subscript and the diffusion coefficient 
for particles can called with appropriate subscript when specific di- 
rections are mentioned. 

Vanoni also found that appears vary with the particle concentration, 
yet, there little information this variation that cannot con- 
sidered the calculation removal. However, must important ef- 
fect when large. hoped that future research will yield quantitative 
information about the effect this paper, the effect ignored. 


Equation Continuity 


The motion particles due settling, mean temporal fluid velocity and 
turbulent diffusion produces flux particles one all the x,y, and 
directions. The resultant flux can represented the vector. order 
describe the vector let and unit vectors the x,y and z-direction 
respectively and let positive the direction positive 

Consider first the flux i-particle. The resultant represented the 
vector, aj, where 


The conservation i-particles x,y,z,t dictates that the divergence vector 
plus the time rate change minus any source flow i-particles 
zero. Hence, 


af. 
where 


The term, Pj, Eq. (3) represents distributed source i-particles 
y,z. This term accounts for the effects hindered settling and flocculation. 
particle class becomes attached another the process floccu- 
lation will generally experience change settling velocity. Hence, dis- 
appears from class and appears another class. Similarly particle may 
have one settling velocity when the local concentration has certain value; 
when the local concentration changes, hindered settling may change the ve- 
locity wj. follows that must represent the rate which particles ac- 
quire settling velocity less the rate which i-particles acquire other ve- 
locities. The units are mass per unit time per unit volume. 

For incompressible fluid the expansion Eq. (3) becomes 
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SUSPENSIONS 


Except for the last term this equation similar that expressed 
McNown,(4) van Driest,(5) and others. 
There Eq. (5) for each value summing these equations and 


introducing Eqs. (1) and (2) into the sum, the continuity equation for total 
concentration becomes 


Since flocculation and hindered settling produce mass, 


(7) 


Substituting Eq. (7) into (6) and assuming that the same for all gives the 
following equation: 


The effects flocculation and hindered settling not appear explicitly 
Eq. effects appear implicitly, however, because they cause 
changes Eq. (8) cannot solved without information how will 
change with particle concentration, flocculation, differential settling, and dif- 
fusion. This information will obtained from studies the settling proper- 
ties the suspension. short, the settling properties the suspension oc- 
cur implicitly Eq. (8) the term 


Initial and Boundary Conditions 


The solution Eq. (8) requires knowledge initial conditions bounda- 
conditions both. becomes necessary, therefore, think about specific 
physical situations. For present purposes, consider the section open chan- 
nel shown Fig. The distance along the channel called and taken 
positive the direction flow. The distance down from the free surface 
the suspension called and positive downward. The y-direction 
normal and such way produce right-handed xyz system. 

For this channel, the initial condition the settling velocity distribution 
point and obtain the settling velocity distribution the sample ex- 
perimental analysis. the distribution not constant over the channel cross- 
section, samples must taken from several points the section and each 
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sample must analyzed. The result will spatial distribution settling 
velocity distribution and this result the initial condition for settling 
the channel. 

The boundary conditions occur the top surface the suspension and 
the bottom the channel. the top, the condition the rate which parti- 
cles cross the surface. Usually particles pass through this surface. But, 
density current along the bottom reservoir settling tank, material 
may settling out relatively still water above into the current below. 

the bottom the channel, the boundary condition depends the behavior 
the particles which have already settled the bottom. particles which 
once settled the bottom the channel are picked into the flow, the 
process called entrainment resuspension. Let the rate which parti- 


cles are picked per unit area bed The boundary condition can then 
stated 


(9) 


This equation the form used describe pickup uniform 
particles turbulence jar. Dobbins verified this equation experimentally 
Theoretical work Lane and shows that will depend upon 
the settling velocities the settled particles, the roughness the bed 
pattern turbulence near the bed. Thus for given flow and bed conditions, 
will depend the properties the suspension. 
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SUSPENSIONS 


The resuspension particles must studied experimentally for each sus- 
pension. During the study, the flow conditions should varied and resus- 
pension measured. The end result should give the resuspension rates 
function flow conditions for the suspension being studied. Experiments 
this sort are currently progress the California Institute Technology. 

appears that great deal research still done before the results 
can used calculating removal. 


Solution the Continuity Equation 


Except for few simple cases necessary solve the equation con- 
tinuity some numerical method. For this purpose, Eq. (8) written the 
form 


Eq. (10), ex, and are all determined measurements. 
particle concentration, the e’s will also function such case the 
relationship between and will have determined experimentally. 

Eq. (10) form easily reduced finite-difference form. For this 
purpose, the space occupied the suspension divided into rectangular 
parallelpipeds shown the upstream end the channel Fig. Each 
parallelepiped has the dimensions Ax, and shown. The corners 
the parallelepipeds form lattice points where the value found. 

The subscript will indicate that point the lattice has x-coordinate 
Ax, while subscripts and indicate coordinate mAy and nAz the 
and z-directions respectively. Using this notation, Eq. (10) can written 
the following form, for steady state settling. 
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use Eq. (11) necessary know and all points the plane 
This knowledge constitutes the initial condition. Any resuspension 
particles between the planes and constitutes the boundary con- 
dition and must also known. With these known quantities, possible 

order repeat the process, necessary know points the 
plane Hence, the fundamental part this calculation determine 
how changes between the plane and or, general, between any 
two planes, This change will depend the settling proper- 
ties the suspension. 

Once the calculated over cross sections various values 
simple matter compute the rate which particles are removed any 
distance, rAx. Thus the problem stated part (a) above can solved. 

The procedure outlined here requires more information than usually 
available the engineer. Its main value, present, showing specifical- 
what research and development are necessary for complete calculations. 
Methods measuring fluid velocity and turbulent diffusion must improved, 
methods measuring resuspension must developed, and methods studying 
the settling properties suspension must expanded. 

Furthermore, the finite-difference calculation shows that information about 
the settling properties suspension must obtained very specific 
form. This information must relate changes the local mean settling velocity 
W(x,y,z) particle settling, particle concentration, flocculation, hinderance, 
turbulence and any other conditions prevailing the immediate neighborhood 
the point x,y,z. program research was initiated for the purpose ob- 
taining this kind information about suspensions. The progress made date 

Due the problems involved calculating removal, natural think 
studying settling hydraulic scale models. However, the model study may 
misleading unless the similarity settling considered. Such similarity 
discussed next. 


Scale Models and Similarity Sedimentation 


The relationship between settling the model and the prototype best ob- 
tained reducing Eq. (8) dimensionless form. All velocities the equation 
can given terms characteristic velocity and dimensionless 
velocities U*, and W*. The result 
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SUSPENSIONS 


Similarly, distance can written terms characteristic length and 
diffusion coefficients can written terms characteristic coefficient, 
follows: 


= e = 
dimensionless concentration can also written terms character- 
istic value, 


Dimensionless time can written terms characteristic length and ve- 
locity. 
Substituting these dimensionless quantities Eq. (8) gives 


There are three kinds similarity involved Eq. (12). The first two are 
geometric and kinemetric. Because them the quantity will the 


same dimensionless constant model and prototype. Furthermore, the value 
the prototype. Hence geometric and kinematic similarity demands that the 
settling velocities the suspension scaled the same ratio the fluid 

The third type similarity concerns the changes settling velocities. 
Suppose that the prototype, the mean settling velocity changes amount 
some distance Ax. For similarity must change the model such 
that the dimensionless rate change the same for both model and proto- 


type. Denoting scale model quantities and prototype quantities 
follows that 


or 
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If, like most open channel models the model was designed according the 
Froude model law, Eq. (13) becomes 


(14) 


The flocculation may related turbulent mixing which scaled according 
Froude’s law follows 


3/2 


Consequently the ratio the effect flocculation equal the negative half 
power the length ratio while the ratio mechanism contributing floccu- 
lation equal the three-halves power the scale ratio. 

Eqs. (14) and (15) show that the response the suspension flow con- 
ditions cannot the same model and prototype. follows that suspension 
with “scaled” settling properties must used the model. Usually, how- 
ever, the same suspension must used model and prototype. this situ- 
ation there are two approaches using models. 

The first approach use the model predict only the flow the proto- 
type; separate experiments are performed study the settling properties 
the suspension. The measured settling properties are then used calculate 
the removal that will occur the predicted flow. This method suited the 
study such problems silt deposition reservoirs, sedimentation the 
mouths rivers, and the diffusion and sedimentation sewage other 
wastes bays and estuaries. 

the second approach the settling studied the scale model. The 
model results are then “scaled up” predict prototype results. order 
perform this scale-up, necessary develop scale equation giving the 
removal the prototype function removal the model. From the 
foregoing analysis obvious that such equation should based much 
the properties the suspension the model laws. Therefore separate 
experiments the settling properties should conducted conjunction with 
the model tests. This second approach frequently used pilot plant studies 
settling tanks. 


Research Quiescient Settling 


Purpose and Scope 


The purpose the research described this section was develop 
methods for determining the settling properties individual suspension. 
More specifically, the objective was measure the settling velocity distri- 
bution suspension, the local mean settling velocity various points the 
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SUSPENSIONS 


suspension, and the factors affecting the local mean velocity. Most the 
work was confined quiescent settling. 

The approach was both experimental and theoretical. The experiments con- 
sisted allowing suspension settle quiescently vertical tube. During 
the settling, small samples were withdrawn from various locations the tube, 
and these were analyzed for suspended-solids concentration. supplement 
the experiments, theoretical analyses were made the quiescent settling 
discrete particles settling tube, the kinetics flocculation during quies- 
cent settling, and the analyses data from the experiments. 


Settling Velocity Distributions 


The purpose the first phase the research was measure the settling 
velocity distribution many varied suspensions. The analysis(7) was 
used for the measurement. portion suspension was shaken and allowed 
settle vertical glass tube. During settling, samples were withdrawn 
from known distance below the top surface the suspension means 
broken-tip pipette small glass tube. 

The results some the experiments are shown Figs. The 
ordinate represents the ratio the particle concentration the sample the 
average concentration the portion the beginning settling, while the ab- 
scissa indicates that the sample was taken depth below the surface time 
after the beginning settling. 

Under the following conditions, curve from these figures represents 
cumulative frequency distribution settling velocities. 


(1) the beginning settling the particles each class settling veloci- 
ties were uniformly distributed throughout the suspension. 

(2) flocculation, turbulence, hindered settling occurred during the 
experiment. this case, the settling velocity distribution the major 
settling property the suspension and often sufficient for calculating 
removal. and have calculated the removal for 
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certain cases particles settling without flocculation hindered sett- 
ling. 


Fig. shows the results pipette analyses performed primary effluent 
from two large sewage treatment plants, primary and secondary effluent from 
third, and various mixtures effluent with digested sludge and sea water 
with digested sludge. The analyses were made the writer connection 
with studies the marine disposal sewage and sludge the vicinity Los 
Angeles, California.(8) These were made apparatus which maintained 
the temperature the suspension 31.5° The average depth sampling 
was about cms for each analysis. 

Even though the pipette analysis takes account flocculation hindered 
settling, the curves Fig. indicate large differences between the settling 
properties the various suspensions. For example, curves and all 
represent primary effluents from sewage treatment plants. However, the ap- 
proximate median settling velocities suspensions and differ factor 
10. 

Pipette analyses were also performed raw Pasadena sewage obtained 
from trunk sewer the Los Angeles County Sanitation Districts. This 


TABLE 


Sources the Suspensions Represented Figure 


Initial Concentration 
Curve Source (equal Solids) 


Primary effluent from the treatment 893 
plants the Orange County 
(California) Sanitation District 


Primary effluent from the Hyperion 212 
Sewage Treatment Plant (City 
Los Angeles). 


Primary effluent from the Joint 314 
Disposal Plant the Los Angeles 
County Sanitation Districts 


Primary effluent plus volume 527 
digested sludge from the Joint 
Disposal Plant, LACSD 


Primary effluent, plus elutriation 289 
effluent from Hyperion Sewage 
Treatment Plant 


One part digested sludge from 
Hyperion Sewage Treatment Plant 
plus parts sea water 


Secondary effluent from Hyperion 
Sewage Treatment Plant 
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sewer serves most Pasadena, San Marino, South Pasadena, and parts 
contiguous communities, with total sewered population about 200,000. The 
samples obtained from the sewer were gross samples; they were taken the 
morning and tested the afternoon the same day. Before the experiment, 
the gross sample was poured into large ceramic crock. When smaller 
portion was required for testing, the gross sample was stirred, and while 
stirring continued, the smaller portion was taken. 

five-gallon gross sample was obtained each several days, anda 
portion was taken from each for pipette analysis. General information about 
the analyses given Table and the results are plotted Fig. These 
curves show variation settling properties the sewage from one run 
the next. For some engineering purposes, would desirable represent 
the data from all four runs single curve. The median settling velocity 
for such curve would about 0.02 cm. per sec., and the median velocities 
for the individual curves differ from this value factor 1.5 less. 

This deviation small compared the differences between median velocities 
Fig. Hence, for the purpose comparing Pasadena sewage with the sus- 
pensions Fig. all four runs could represented single curve. 

From the curves Figs. and two conclusions can drawn. First, 
suspension which may considered similar can have significantly 
ent settling properties. Second, even for suspension heterogeneous 
sewage, the difference between runs for single suspension can small 
compared the difference between suspensions. 


Table 


Pipette Analyses Pasadena Sewage 


Settling tube liter graduate 
Volume used liter 
Temperature control none’ 

Size samples withdrawn ml, 


Method withdrawing samples ml, broken tip pipette lowered 
into suspension hand for each sample, 


Depth samples below surface sewage, 


Gross Pipette Initial Temp. 
Sample Analysis Concentra- During 
Date Obtained Began tion Test 
Thur. 21/54 8:30 2:35 304 
Mon, 15/54 9:00 1:45 374 


Mon, Nov, 22/54 9:00 1:00 380 23-30°C 
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Part the difference between the curves Fig. may due sampling 
error. Whenever portion withdrawn from heterogeneous suspension 
such raw sewage, the properties the portion may not the average 
properties the suspension. test this error, tests were made two 
portions from single gross sample. means two identical glass settling 
tubes placed single constant-temperature water bath, two similar portions 
sewage from single gross sample were subjected identical pipette 
analyses. The portions were approximately four liters volume. 

The experiment was rerun for second gross sample, and the results 
both runs were plotted the graph shown Fig. The difference between 
curves for one run nearly large the difference between curves shown 
Fig. Hence, some the variation the latter curves due sampling 
errors. 

Sampling errors can reduced increasing the volume the portion 
used the suspension. The portions used these pipette analyses varied be- 
tween and liters, and they appear too small. Hence larger samples 
are recommended for suspensions considerable heterogeneity. 


Factors Affecting 


Unless the particles suspension settle without flocculation hindered 
settling, the pipette analysis gives only rough indication the settling 
properties suspension. more general situation necessary use 
experiment that will show how flocculation and hindered settling affect the 
settling velocities the particles. 

hindered settling the concentration suspended particles high enough 
for the presence one particle affect the settling its neighbors. 
flocculation, the other hand, faster particle overtakes slower one and 
becomes attached it. The two settle henceforth unit with velocity 
usually different from the original velocity either particle. These two ef- 
fects cannot separated experiment. they occur simultaneously the 
result single change However, order find out precisely how 


c 

= 

2a 


Pipette analyses Pasodeno sewage(See Toble2). 
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they produce the changes convenient think them separately. 

The hindered settling uniform particles without flocculation has been 
studied theoretically and experimentally McNown and Lin(9) and 
Steinhour, (10) among others. The former found that the individual settling 
velocity particle depends upon the volume concentration the particles 
and upon the Reynolds number the particle. Steinhour, the other hand, 
produced formula giving velocity function volume concentration alone. 
This formula was based experiments where particles were settling the 
stokes range. 

Both these investigators considered cases where the particle concen- 
tration was uniform throughout the settling tube. Kynch(11 went further and 
made theoretical analysis uniform particles tube where the concen- 
tration varied with depth. assuming that the individual particle velocity 
depended solely the concentration the neighborhood the particle, 
obtained results which agree with experimental observations the subsidence 
thick slurries. 

the basis these studies was concluded that the ratio the settling 
particle velocity secondary factor. Hence, suspension with particles 
many velocities, change will cause the velocity each particle 
change the same ratio. However, the magnitude the change velocity 
will equal the product the ratio and the particle velocity. Thus, the magni- 
tude will increase with particle velocity. 

follows, that for given change the mean settling velocity will 
change amount depending upon the distribution settling velocities 
the time the change. This distribution can probably characterized the 
mean settling velocity and the standard deviation, the velocities. There- 
fore, changes will function itself, and the volume concen- 
tration particles. 

For flocculation theories existed which were suitable for this research. 
The previous theories were all confined predicting the rate interparticle 
contacts, while ignoring the effect settling velocities. Furthermore, they 
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Tube 


Initial conc. 436 326 


emperature 29.5 295 


0.01 
Fig. Duplicate pipette analysis Pasadena sewage. 
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dealt with suspensions having velocities only few classes. Therefore, the 
kinetics flocculation during quiescent settling was studied theoretically(12) 
for the purpose discovering how flocculation produces changes Be- 
cause the number unknowns the theory could not completed, but was 
carried far enough show that the change point depends the 
following properties the suspension the same point. 


(1) The volume concentration particles. 

(2) The mean settling velocity particles. 

(3) The standard deviation the settling velocities the particles. 

(4) The fraction inter-particle contacts that result the union the 
particles. 


Item (4) simply means accounting for the surface chemistry the 
particles. 

The value point will also change free settling. That is, the 
absence flocculation and hindered settling easily shown(12) that 


Thus, the change depends itself and the standard deviation 
settling velocities. 

appears, therefore, that the first three the items listed above con- 
nection with flocculation are also the most important items connection with 
free and hindered settling. Consequently, experimental study quiescent 
settling should involve the determination least and concentration 
volume. The determination the first these, described part (d) 
which follows. 


Measurement 


Since measurement has not been reported the literature, some 
time was spent devising method. Finally, was decided use ex- 
periment based suggestion The suspension was allowed 
settle quiescently vertical tube pipette analysis. During the sett- 
ling series samples were taken each several depths, and these sam- 
ples were analyzed for particle concentration. This experiment was given the 
name multiple-depth pipette analysis or, simply multiple-depth analysis. 

Some very simple equations will show how such analysis can used 
measure 

For the purpose discussion, consider the hypothetical settling tube shown 
the left Fig. All horizontal cross-sections the tube are unit area, 
and any time, the particle concentration, over any such cross-section 
constant. Hence, function the depth below the top the suspension, 
and the time after beginning settling. 

the beginning the test pipette samples are withdrawn simul- 
taneously many depths. The results are plotted the curve, the 
and the results are plotted the curve, The curves, andt 
represent similar operations times and respectively. 

Physically speaking, each these curves profile the concentration 
some given instant. When all the profiles from analysis are put one 
diagram, can called the concentration profile diagram. Mathematically 
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The concentration profile diagram used conjunction with the continuity 


SUSPENSIONS 


ter. any stage the discussion, possible give the diagram physi- 
cal mathematical interpretation, depending which more illuminating. 


equation for the calculation W(z,t). First, noted that for the settling 
tube, Eq. (8) reduced the following: 


Integrating this equation with respect gives 


particle concentration. 


z=depth below surface. 


Settling Concentration profile 
tube diagram 


Fig. Analysis data from multiple-depth 
pipette analysis. 
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Eq. (18) shows how can calculated. value selected; for ex- 
ample, consider Fig. The areas 025, 024, and 023 are calculated and 
the values are plotted against this done for sufficient number pro- 
files, the result will smooth curve giving the area under the profile 
above function time. The slope the curve precisely the right 
hand side Eq. (18), and hence equal Dj. The process can 
repeated for any value 
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zero forall depth 
measurements 
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surface reading 
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Glass settling tube 


Sampling No. 

intake 

(Glass 
4mm. 


Fig. Multiple- depth settling tube. 
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order test the feasibility this approach measuring pilot ex- 

was performed using the apparatus shown Fig. Sampling in- 
takes, pipettes, were located three depths shown. Samples could 

drawn through the top two means vacuum and through the bottom one 
means gravity flow. 

The suspension tested was thoroughly mixed separate container 
and then poured quickly into the tube. Since the upper two intakes are attached 
platform which rests freely the tube, these intakes could put 
place immediately after pouring. After they were place, the elevation the 
top surface was read the scale shown and the withdrawing samples began. 

The surface elevation was read again after the withdrawal each sample. 
With these readings was possible allow for the lowering the surface 
estimating the distance settled particles. The estimated distance was used 
for each sample. Time, was measured from the end pouring. 


Table 
Multiple- Depth Analysis Suspension 


Clay and Alum in Water 


Time Depth Weight 
After of Volume Temp. of % of 
Sample Sample Sample Conc, 
180 3 103, 2 0, 57 33 21.6 104 
; 245 2 39. 7 0. 16 28 18, 1 103 
300 l 18, 4 0.062 37 23.6 102 
360 3 101.4 0, 28 31 1/2 20, 0 101 
| 480 1 i? 1 0, 036 28 16,7 95 
600 2 37.5 0, 063 331/22. 30.3 21.4 108 
j 720 3 99.2 0. 14 32 30, 3 22.9 14 
960 1 15,2 0.016 22 8.1 59 
1200 2 35.6 0, 030 25 11,2 71 
1920 1 13, 2 0, 069 27 7.7 45 
2400 2 33, 7 0,014 31 9.0 46 
2880 3 95. 4 0.033 29 31.5 s Ae 40 
3840 ll, 0, 0029 31 32.0 26 
5760 3 93. 3 0.016 27 32,8 4.6 27 
7740 1 8.9 0.00i2 31 33.4 2.7 14 
10100 2 29. 4 0, 0029 40 34.0 4.1 16 
14100 3 91.1 0, 0065 26 34.8 1,2 7 
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This apparatus was used for the multiple depth analysis suspension 
bentonite clay Pasadena tap water and alum. The 
initial concentration clay was 655 mg. per liter, while that alum was 
mg. per liter. The results the experiment are given Table and plotted 
the concentration profile diagram Fig. Since samples were taken 
only three depths some intermediate plotting and interpolation were necessary 
produce the profiles. 

Once the profiles were obtained, the mean settling velocity W(z,t) was 
calculated the manner described above. Values for three values are 
given Fig. possible calculate the order magnitude 


(z,t)=% weight initial concentration 


100 
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depth below surface incm. 
—, 
5000 


100 
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Intake 
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Initial concentration=655 


Fig. Concentration profile diagram, Bentonite 
clay and alum water. 
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flocculation effects local settling. depth cm. the mean settling 
velocity increases from 0.021 cm. per second 200 seconds 0.037 

500 seconds. This change represents increase per cent and 
average rate increase 0.000053 cm. per sec. per sec. Furthermore, the 
change concentration during the same time was only per cent. The effect 
particles settling out suspension must have been small. Therefore, the 
rate change primarily due flocculation. 

possible tell directly that flocculation causing the particles 
speed looking the dashed lines constant z/t Fig. The physical 
significance these lines best explained the following manner. ob- 
server the surface the suspension and descends through 
the suspension constant velocity. The concentration that observes 
various depths given line z/t equal the velocity. 

can shown (12) that when neither hindered settling nor flocculation oc- 
cur, these lines are straight and parallel the z-axis. hindered settling 
slows the particles down more than flocculation speeds them up, the lines 
slope away from the z-axis depth increases. the converse true, the 
lines slope toward the z-axis depth increases. For Fig. then, seen 
that flocculation the predominant effect. 

The maximum mean settling velocity for cm. occurs 700 seconds. 
Its value 0.056 cm. per second 2.8 times the mean settling velocity 
200 seconds. The average rate change between 200 seconds and 
700 seconds was 0.00007 cm. per sec. per sec. 

maximum value occurs all depths the settling tube. After this 
value, the mean settling velocity decreases. Particles are settling out sus- 
pension, and the loss faster particles offsets the effect flocculation. 

While this work was being done, Fitch(14) presented the results 
water. The concentration was 400 parts per 
million ppm), while the concentration ferrisul was ppm. This suspension 


mean settling velocity cm. per sec. 


t=time sec. 
Bentonite clay and water. 
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was allowed settle tube seven feet deep and 5-1/2 inches internal 
diameter. During the settling, samples were taken seven depths means 
veterinary hypodermic needles which passed through the walis the tube. 
For temperature control, the outside the tube was covered with insulation 
one inch thick. 

Fitch’s results are presented concentration profile diagram Fig. 
The solid curves are concentration profiles lines constant time and 


depth below surface feet 


Z= 


Fig. Concentration profile didgram, and 


water 
(Based data Fitch) 
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the dashed curves are curves constant z/t. For >0.25 hours and 
feet, the constant z/t lines slope toward the z-axis, showing the effect 
flocculation. 

Following the method described above, the local mean settling velocity was 
calculated. The mean velocities for feet and feet are shown 
Fig. 10. feet, the mean settling velocity increased from 0.013 cm. 
per sec. 0.25 hours 0.035 cm. per sec. 0.5 hours. The average 
rate change during this time was 0.000025 cm. per sec. per sec. 


hours, the concentration feet per cent the initial concentration. 


Consequently, this time the change caused primarily floccu- 
lation. 

interesting compare the change Wfor whiting and ferrisul with 
the change for clay and alum. The comparison valid only for those stages 
settling for which the concentration has not decreased significantly. For 
whiting and ferrisul, the rate change during this stage settling was 
0.000025 cm. per sec. per sec. For clay and alum, the corresponding rate 
change was 0.000053 cm. per sec. per sec. The values differ factor 
two. 

One important practical conclusion can drawn from these determinations 
That that the effect flocculation increases with depth from the 
surface. This shown the curves for Figs. and 10. the settling 
tube either case had been deeper, the peak values would probably have 
increased with depth until limiting value was reached. 
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Fig. 10. Local mean settling velocity function 
time (based data 
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Proposals for Future Work 


The next logical step this work design experiment which 
possible measure and particle concentration. For such experi- 
ment, desirable have tube about ten feet deep with sampling 
intakes. The sampling intakes should enter through the walls the tube 
cause minimum obstruction settling particles. The tube should 
housed water jacket for temperature control, but both jacket and tube 
should have windows for visual observations. 

The cross-sectional area this tube should least equal that 
circle twelve inches diameter. Because this large area, the samples 
withdrawn during settling can large 500 ml. Then each sample can 
allowed settle 500 ml. graduate for simple pipette analysis. 

Five six ml. samples from each 500 ml. graduate should enough 
give settling velocity distribution for each 500 ml. sample. From this 
distribution possible calculate the mean and standard deviation the 
settling velocities for each sample. course, flocculation may occur the 
graduate, but the effect small compared that the tube. 

The same 500 ml. samples can used for determination and the 
concentration volume. From the former the concentration profile can 
drawn. The values can then calculated shown above and checked 
against the value obtained from the simple pipette analysis. 

Before this more elaborate multiple depth-analysis could performed, 
the experiments settling analysis were discontinued, and, yet, pro- 
vision has been made for their continuance. hoped that someone will find 
the problem sufficient interest continue the work. so, will eco- 
nomical modify the tube allow for the study turbulence. 

matter extend the proposed multiple-depth analysis in- 
clude effects turbulence. Rouse(15) and Dobbins(3) have shown that uni- 
form field turbulence can created placing vibrating grid the 
settling tube. Furthermore, Dobbins devised method for controlling the rate 
particle resuspension the bottom the tube. these additions are 
corporated into the apparatus proposed above, will possible study 
flocculation turbulent fluid. 


The Multiple-Depth Analysis Approximate Calculations Removal 


independent depth. Then, the settling particles similar settling 
vertical tube which moving velocity the x-direction. The pro- 
files the tube time, corresponds profiles the channel distance 
follows that these profiles can used directly calculate removal. 


Short-Cut Methods Calculating Removal 


“Ideal” Settling Tanks 


length width and depth making certain simplifying assumptions 
about such basin, devised the concept “ideal tank”. Using 


Part (a) Fig. shows the settling zone rectangular settling tank 
the nomenclature Section these assumptions can stated follows: 


U(x,y,z) constant 
V(x,y,z) W(x,y,z) 
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resuspension the bottom the tank 


Although these assumptions are sufficient describe the tank, Camp also 
assumed that all particles the tank settle without flocculation hindrance. 


Since this last assumption unnecessary, this paper the term ideal tank 


will refer only the assumptions about the tank. place the assumption 


“<i 


(a) Rectangular 


(b) Radial flow 


Settling zones ideal tanks. 
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about the settling, the results multiple-depth analysis can used. This 
approach based the fact that ideal tank settling identical that 
vertical tube, and concentration profiles distance the tank corre- 
spond profiles time the tube. the concentration profiles from the 
tube are used, the removal the ideal tank will given 

where 
and ¢(z,t) represents concentration the tube. 
settling tanks, the removal ratio, r(x), more important than the re- 
moval itself. The former defined the fraction particles entering the 
tank that settle out distance Hence, 
Oo fe) 
The integrals Eqs. (19) and (20) can represented the concentration 


profile diagram areas. Consider, for example, the diagram shown Fig. 
For tank depth and distance the removal represented 
area 0561 while the removal ratio represented 


area 0561 
similar approach can used for radial flow tank. Part (b) Fig. 
shows what might considered “ideal radial flow tank”. Fluid enters the 
tank through the surface cylinder radius centered along the vertical 
axis the tank. then flows radially and horizontally the cylindrical walls 
the tank. These walls have radius 
The flow conditions for this tank are assumed 
x y Z | d 
The initial conditions are assumed 
while the bottom boundary condition the same that for Camp’s 
Settling this radial tank identical settling vertical tube. 
Concentration profiles radius the tank correspond profiles time 
(21) 
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the tube. follows that removal the tank radial distance, given 


R(P) 


where given Eq. (21). Eq. (22) shows that removal the radial tank 
represented area the concentration profile diagram just for 
the rectangular tank. 

suggested that the combination assuming ideal tank and perform- 
ing multiple-depth analysis the suspension will give first approximation 
the behavior the tank. This approximation will often better than mak- 


ing elaborate study the hydraulics the tank while ignoring the proper- 
ties the suspension. 


Ideal Tanks with Turbulence 


For better approximation rectangular cylindrical basin, 


possible add turbulence the ideal tank. the turbulent diffusion as- 
sumed 


e 


constant, 


the settling the tank will equivalent settling tube into which uni- 
form turbulence introduced. 

and have assumed this sort tank study the 
effect turbulence retarding settling. Camp assumed that the particles 
settled without flocculation hindrance. However, this assumption un- 
necessary, simple matter put the actual suspension settling 
tube deep the tank and introduce turbulence means vibrating grid. 

The combination assuming ideal tank and performing turbulent 
multiple-depth analysis the suspension suggested second approxi- 
mation the behavior rectangular basin. For better approximations, 


necessary obtain enough information for the step calculation outlined 
Section 


Overflow Rate and Detention Time Ideal Tanks 


the technical literature, one finds great deal discussion about 
whether settling tank should designed the basis overflow rate de- 
tention time. The discussion about the merits these two approaches in- 
conclusive, because disregards the properties the suspension. will 
now demonstrated that whether overflow rate detention time should 
used depends upon the nature the suspension. Indeed, for single tank, the 
removal ratio may closely related overflow rate for one suspension, 
detention time for second, and neither for third. 


The detention simply the average time that element fluid remains 
the tank. Let called For the rectangular ideal tank, 


(23) 
while for the radial, 
(24) 


Overflow rate, the other hand, obtained dividing the volume rate 
flow through the tank the horizontal area the tank. Let For the 
rectangular basin, 


(25) 


while for the radial, 


(26) 


Eqs. (25) and (26) show that for constant flow rate and overflow rate, de- 
tention time varies with depth. Thus, the effectiveness the tank de- 
pend overflow rate, r(x) will not vary with long D/T constant. 
Conversely, the effectiveness depends detention time, r(x) will not vary 

study the problem experimentally, multiple-depth analysis per- 
formed the suspension concerned. The settling tube for this analysis 
should deep the deepest possible tank considered. The results 
from the analysis are plotted concentration profile diagram shown 
Fig. 

this figure, let and represent the depths three tanks 
considered. The detention times for these tanks are and T3, re- 
spectively. The tanks all have the same overflow rate. Therefore, 


w = 


and the dashed line line constant z/t. The other dashed lines 
are also lines constant z/t. 
the removal ratio depend overflow rate alone, then 


area 0561 area 13, 16, (27) 
area 0261 area 11,1 area 12, 


Furthermore, this type relationship must hold for any overflow rate, i.e., 
for any dashed curve constant z/t. can that Eq. (27) will 
hold for arbitrary overflow rate and only the lines constant z/t are 
straight and parallel the z-axis. 

If, the other hand, removal ratio depends only detention time, 
equation the form 
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must hold for each profile. Eq. (28) will hold and only the profiles are 
all straight lines parallel the z-axis. 

follows that for quiescent settling ideal tanks the problem overflow 
rate and detention time determined the pattern the profile diagram. 
This pattern, turn, determined the settling properties the sus- 
pension. 

example suspension for which detention time more important, 
consider the data obtained Fitch and plotted Fig. above. The profiles 
are almost straight and vertical. Thus, for any tank six feet depth, 
the removal affected more detention time than overflow rate. Calcu- 
lations Fitch substantiate this conclusion. 

With only these curves hand, one might ask the following question. What 
will the situation when this suspension settles rectangular tank ten 
feet deep? answer this question, necessary try sketch the 
lines constant z/t between the depths six and ten feet and draw the pro- 
files from these. After few trials, becomes evident that the profiles are 
not likely change their shape drastically. Hence, removal will still depend 
primarily detention time. 


Detention Time and Overflow Ratio—General 


The data Fitch shows that when flocculation occurs ideal tank the 
removal ratio can more closely related detention time. Furthermore, 
reference (12) shown that flocculation has tendency cause the verti- 
cal straight profiles which indicate the dependence detention time. 

The reason quite simple. When particles settle without flocculation, the 
concentration given time will normally increase with depth. Hence, the 
profile slopes away from the z-axis. With flocculation, however, the faster 
particles settling through the slower ones gather the slower ones and drag 
them out suspension. This tends cause the concentration any given 
time decrease with depth. The combined result can cause the concentration 
independent depth, i.e., straight profile parallel the z-axis. 

However, has shown that when the particles settle without floccu- 
lation hindered settling, removal depends overflow rate alone. Between 
the suspensions discussed Camp and those discussed Fitch are sus- 
pensions for which both overflow rate and detention time are significant. The 
diagram Fig. represents example this last situation. 

When turbulence during settling considered, the studies overflow 
rate and detention time should based results from turbulent settling 
analysis. yet, the available experimental results this type are too mea- 
ger permit any conclusions. 


Sediment Deposition Canals 


the design canals, necessary consider minimum velocities for 
non-silting. the latter cannot maintained, desirable know where 
the canal the silt will deposited. Some useful information about this 
problem can obtained from turbulent multiple-depth analysis. 

This analysis should performed tube deep the canal. Inside 
the tube, the level and vertical distribution turbulence should close 
possible that the canal. The vertical distribution turbulence can 
obtained varying the grid spacing and bar size along the vibrating grid. 


The turbulence level can obtained adjusting the frequency the 
vibrations. 
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The tube can used find out what velocity necessary for keeping the 


particles suspension. letting the suspension settle the tube with vari- 


ous grid frequencies, will possible decide upon the turbulence level 


that keeps deposition reasonably small amount. The velocity the canal 


which produces the turbulence the necessary velocity the canal. 

For each grid frequency, the test should run for considerable time. 
After only short time the tube, the particles suspension may 
reasonably well suspended given level turbulence. However, floccu- 
lation occurs, the particles may increase the settling velocities and begin 
settle out. This process may not evident for thirty minutes hour. 

minimum velocity cannot maintained the canal, the results the 
analysis will show the pattern silt deposition. For the flow the canal 
fairly uniform, profiles distance the canal correspond profiles 
the tube. Removal calculated the manner described for the ideal tank. 

When the level turbulence the canal sufficiently low, the results 
quiescent multiple-depth analysis may give some information about depo- 


sition. Fig. represents the results quiescent analysis bentonite clay 


depth below surface feet. 


Calculated from pipette analysis 


Measured multiple-depth tube. 


clay. 


Initial concentration= 
Temperature= 26°C 


Fig. Concentration profiles for bentonite 
clay Pasadena tap water. 
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SUSPENSIONS 


Pasadena tap water. The analysis was performed lucite tube with 
internal diameter 3-3/4 inches and depth feet. Thus the settling 
might correspond, roughly, settling during tranquil flow irrigation 
canal feet deep. 

preparing the suspension, the clay was blended thoroughly with 0.5 liters 
water. This blended mixture was then shaken with enough water bring 
the final clay concentration 872 per liter. After two minutes vigorous 
shaking the suspension was poured into the tube. This mixing might simi- 
lar that which the suspension would receive passing through diversions 
works into the canal. 

During the settling, samples were taken four depths. The concentrations 
these samples were used plot the profiles shown solid lines Fig. 12. 
For the purpose comparison, the settling velocity distribution the origi- 
nal clay particles was obtained means pipette analysis suspension 
clay and deflocculating agent distilled water. From the settling distri- 
bution, the dotted profiles Fig. were determined. These profile show 
how the particles would settle without flocculation. 

the basis the dotted profiles, one would expect that the clay would 
carried suspension indefinitely. Therefore would pass through the canal. 


fact, however, most the clay would have settled the bottom the canal 
one hour. 


CONCLUSIONS 


The principal conclusions based the research are follows: 


(1) order predict the settling particles flowing suspension, 
necessary know the properties the flow and the settling properties 
the individual suspension. general, best determine both 
direct measurement. choice must made, will often better 
determine the properties the flow calculation reasonable as- 
sumption, and measure the properties the suspension. 

(2) the initial stages studying suspension, profitable obtain 
settling velocity distribution means pipette analysis compara- 
ble experiment. When flocculation and hindered settling are negligible 
this distribution can used estimate removal; when they are 


not, the results the are still useful comparing suspensions. 


(3) multiple-depth pipette analysis can used study the effect hin- 
drance, flocculation, and turbulence the settling. For quiescent sett- 
ling, the analysis should include measurements the mean settling ve- 
locity, the standard deviation settling velocities and the particle 
concentration various depths and times. For turbulent settling the 
analysis should include these three measurements plus any additional 
measurements related the turbulence. 

(4) Whenever the settling the flowing suspension can approximated 
settling quiescent turbulent settling tube, the concentration pro- 
files from multiple-depth analysis can used calculate removal. 

(5) When suspension flows through settling tank, the removal parti- 
cles may depend upon detention time, overflow rate, both. Whether 

one the other more important depends primarily upon the settling 

properties the suspension. For quiescent settling ideal tanks, 
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multiple-depth settling analysis will indicate the relative importance 
overflow rate and detention time. 

(6) Settling hydraulic model will not similar settling the proto- 
type unless suspension with scaled settling properties used the 
model. When this not possible, two approaches are available. one, 
the model used predict the properties the flow, and separate ex- 
periments are used measure the properties the suspensions. 
These properties are used calculate what the removal will the 
predicted prototype flow. the other approach, the prototype sus- 
pension settles the model. The settling properties the suspension 
are studied separately and these properties are considered using the 
model results predict prototype results. 


This paper based thesis(12) submitted the California Institute 


Technology 1958 partial fulfillment the requirements for the Ph.D. 
degree. The writer grateful Professor Jack McKee for his super- 
vision and generous assistance during the course the research. 


The experimental work was carried out with the aid research grant 


(RG-4405) from the Division Research Grants, National Institutes Health, 
Public Health Service. 


REFERENCES 


Vanoni, Vito A., “Transportation Suspended Sediment Water,” Trans. 
Am. Soc. Civ. Engrs., vol. 111, 1946, pp. 67-102. 


Ismail, Hassan M., “Turbulent Transfer Mechanism and Suspended Sedi- 
ment Closed Channels,” Trans. Am. Soc. Civ. Engrs., vol. 117, 1952, pp. 


409-434. 


Dobbins, William E., “Effect Turbulence Sedimentation,” Trans. Am. 
Soc. Civ. Engrs., vol. 103, 1944, pp. 629-678. 

McNown, John S., Discussion Turbulence Sedimentation,” 
Trans. Am. Soc. Civ. Engrs., vol. 103, 1944, pp. 657-660. 

Van Driest, R., Discussion “Effect Turbulence Sedimentation,” 
Trans. Am. Soc. Civ. Engrs. vol. 103, 1944, pp. 674-676. 

Lane, and Kalinske, A., “The Relation Suspended Bed 
terial Rivers,” Trans. U., vol. 20, 1939, pp. 637-641. 


Krumbein, and Pettijohn, J., Manual Sedimentary Petrography, 
Appleton-Century, New York, 1938. 


Brooks, Normal H., “Predictions Sedimentation and Dilution Digested 


Sludge Santa Monica Bay,” Report Hyperion Engineers, August 
1956. 


Brooks, Norman H., “Settling Analysis Sewage Effluents,” Memorandum 


Hyperion Engineers, July 1956. See also 


Hyperion Engineers, Ocean Outfall Design, Los Angeles, California, 
October 1957. 


j 
4 
| 
4 
a 


10. 


11. 


12. 


13. 


14. 


15. 


16. 


SUSPENSIONS 
McNown, and Lin, P., “Sediment Concentration and Fall Velocity,” 


Proc. 2nd Midwestern Conference Fluid Mechanics, Ohio State Univ. 


1952, pp. 401-411. 


Steinour, H., “Rate Sedimentation,” Industrial and Engineering 
Chemistry, vol. 36, no. 1944, pp. 618-624. 


Kynch, J., Theory Sedimentation,” Trans. Faraday Soc., vol. 48, 
1952, pp. 166-176. 


McLaughlin, Ronald T., “On the Mechanics Sedimentation Artificial 
Basins,” Ph.D. Thesis, Calif. Inst. Tech., Pasadena, California, 1958. 
Camp, Thomas R., “Sedimentation and the Design Settling Tanks,” 
Trans. Am. Soc. Civ. Engrs., vol. 1946, pp. 895-958. 

Fitch, B.,*The Significance Detention Sedimentation,” Sewage and 
Industrial Wastes, vol. 29, no. 10, October 1957, pp. 1123-1133. 

Rouse, Hunter, “Experiments the Mechanics Sediment Suspension,” 
Proc. 5th International Cong. for Applied Mechanics, 1938, pp. 550-554. 


Camp, Thomas R., Discussion Turbulence Sediment- 
ation,” Trans. Am. Soc. Civ. Engrs., vol. 103, 1944, pp. 660-667. 


one, 
{ 
the 
i | 
ealth, 
‘rans. 
edi- 
52, 
Am. 
tion, 
igested 
7 4 
randum 


a 
4 
3 
—_— 
| 
‘ a 
| 
| 
| 
id 
q 
~ 
| 
3 o! 
re 
4 Ir 
4 
q 
a 
= 


2315 December, 1959 


Journal the 


HYDRAULICS DIVISION 


Proceedings the American Society Civil Engineers 


SUGGESTED LEGISLATION FLOOD PLAIN 


Joseph Perrey,! ASCE 


ABSTRACT 


Adequate statutory authority for city, county and regional planning lack 
ing many areas. Such authority must provided communities are 
regulate the use flood plains manner that will limit and tend reduce 
the damage from floods. This problem concerns all levels government. 


SYNOPSIS 


The rapid growth urban and suburban areas has caused rapid expansion 
developments into the flood plains this country. The increase flood 
damageable improvements has outstripped all efforts provide flood control 
flood protection. 

increasing awareness the necessity for limiting the flood damage po- 
tential needed the local level. The preparation plans for the wise use 
flood plains, and the adoption regulations ensure such use are re- 
quired for proper and wholesome community development. 

Not all states provide adequate statutory authority for city, county and 
regional planning. Such authority should provided. addition, states 
should provide controls over construction dams, levees and encroachments 
floodways, should collect basic flood data, and should advise local planning 
agencies regarding flood plain regulation problems. 

Federal agencies should expand basic data collection programs, observe 
local flood plain regulations, and require enactment appropriate flood plain 
regulations local cooperating parties prerequisite qualifying for 
certain types federal programs. 


Note: Discussion open until May 1960. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. Paper 2315 
part the copyrighted Journal the Hydraulics Division, Proceedings the 
American Society Civil Engineers, Vol. 85, No. 12, December, 1959. 

Presented the July 1959 Hydraulics Division Conference Fort Collins, 
Col. 


Chief Engineer, Indiana Flood Control and Water Resources Commission, 
Indianapolis, Indiana. 
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INTRODUCTION 


The rapid growth urban and suburban areas since World War has been 
outstanding feature the economic boom that has taken place during the 
post war years. This growth has been stimulated result the subnormal 
amount building during the depression years the thirties, the war re- 
strictions the early forties, the very rapid increase population following 
the war, and the favorable employment and earning levels that have existed 
during the past fourteen years. 

The land requirements accommodate this growth have been greatly en- 
larged the trend toward sprawling one story houses, and septic tank 
sanitary system drainage requirements suburban locations. These trends 
lead lots two three times larger than those that were considered adequate 
forty fifty years ago. The greater land requirements have placed pre- 
mium areas flat gently sloping topography which can developed 
most economically. Undeveloped flood plains fulfill these requirements, par- 
ticularly the regions the country where topography tends rolling 
hilly. 

For these reasons thousands new homes have been constructed areas 
subject flooding. many instances subdivisions have been developed with 
full knowledge the dangers involved, and the developer has either ignored 
these dangers considered them remote that little consideration needed 
given them. addition, the planning commissions other bodies 
with responsibility for approving new developments frequently have paid little 
attention flood danger element considered appraising the 
suitability site for residential, business industrial development. 
other instances the approving bodies have recognized the danger but have not 
known how deal with it, except bear with and hope for the best. 

Industry bit more cautious than many other developers selecting 
sites for new developments expanding old installations and will try 
avoid flood dangers will protect against them. However, there are numer- 
ous instances throughout the country where such dangers have not been 
avoided. Many times industries have exposed themselves possible damage 
because other considerations that were deemed more pressing and 
significant their operations. 

addition being exposed flood danger, many developments also en- 
croach the floodway, reduce the efficiency and capacity the floodway and 
consequently increase flood levels and damage themselves and others. 

some areas, flood heights have been reduced flood control reser- 
voirs and the danger flooding diminished, developments have encroached 
farther into the flood plains and again have become subject damage from 
lesser floods. Such use the flood plains tends nullify many the benefits 
from flood control projects. 

has been pointed out that although the Federal government 
has spent over billion for flood protection since 1936 there has been 
decrease the mean annual flood loss. the contrary, flood damages have 
tended increase significantly. attributes much this increase the 
great increase use the flood plains. 

the flood damage potential growing and outstripping the relief being 
obtained from flood control and protection projects, why are constructive 
measures not being taken more effectively control the improvements that 
create and increase the damage potential? The answer the lack 
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FLOOD PLAIN REGULATION 


recognition the problem the local level, lack adequate legislative au- 
thority deal with the problem, failure use existing authority, and lack 
adequate knowledge for dealing with the technical phases the problem. 


Local Regulations 


the problem reducing the flood damage potential our flood plains 
considered oneof controlling the unwise actions people, its solution 
will involve the use the police powers the states and the extension 
those powers local governing officials. The use police powers still 
realm that reserved for the states and local governments under the terms 
the Constitution. therefore devolves upon local officials, who are most 
direct contact with the people involved, perform the primary policing func- 
tions. 

There has been growing awareness the responsibility local units 
government plan for the wise use flood plains. During the past five years 
increasing number communities have included flood plains special 
use areas their community plans. However, the number such communi- 
ties pitifully small compared the total number communities engaged 
planning. 

The lack adequate statutory authority permit proper planning 
factor that retarding progress some areas. The use improper tech- 
niques ineffective regulations provides inadequate controls other areas. 

The legislatures states (including Alaska and Hawaii) and the District 
Columbia have expressly authorized community planning. Community-type 
planning regional basis has been authorized states. 

Although planning enabling statues have been enacted states, 
states there express statutory authority granted extend planning 
activities beyond the corporate limits cities; states there ex- 
press statutory authority granted for final approval adoption compre- 
hensive plan the city council other local legislative body after has 
been prepared. 

Some the five states (Arizona, Florida, Missouri, Texas, and Wyoming) 
that not have general laws authorizing the establishment local planning 
agencies with power prepare comprehensive community plans for urban 
areas, have provisions made for planning other means. Florida, author- 
ity establish city and county planning commissions found special 
legislative enactments. Missouri provides for county planning under certain 
circumstances; and Texas, city planning commissions are authorized 
the charters individual cities. 

addition deficiencies planning legislation, there are also deficiencies 
authority for zoning and regulating subdivisions, particularly county 
regional basis. State legislatures have passed enabling acts giving zoning 
powers about half the 3,000 counties the United States.(3) some 
states all counties may zone; others, only the more populous counties may 
so. 

readily seen that the lack adequate authority plan and zone will 
restrict the use flood plain regulations large sections the country. 

the communities that have authority plan and zone, not all are taking 
advantage the powers available them and some are even hostile their 
use. Such hostility usually justified the theory that community planning, 
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zoning and other regulations infringe the right the individual exercise 
free choice the use his property and therefore should opposed for the 
preservation individual liberty. Fortunately, such attitudes are held 
decreasing minority becomes more evident that increasing population 
density requires closer living associations with neighbors and that unrestrain- 
individual liberty frequently produces conflicts between individuals in- 
dividuals and groups people the detriment the general welfare. 

The need for authority plan land use and the need for courage and de- 
termination taking steps guide such use along wise paths indicate that 


there are four primary requirements for carrying out regulation flood 
plains the local level: 


State enabling legislation empower communities prepare com- 
munity development plans. 


State enabling legislation empower communities enact zoning and 
other regulations with flood plain provisions. 


The will prepare the necessary plans and enact the necessary or- 
dinances. 


The will enforce the ordinances after enactment and resist grant- 
ing exceptions them. 


These basic principles appear simple enough but their application usually 


difficult. zoning flood plain serve the needs the community the 


best manner possible, with minimum danger from flood damage and with 
the future growth the community considered, careful analysis the 
flooding situation required. distinction between the general “flood plain” 
and the floodway must made. The flood plain the area adjoining the 
river stream which has been the future may covered flood 
water. The floodway includes the channel river stream and those por- 
tions the flood plains adjoining the channel which are reasonably required 
efficiently carry and discharge the flood water flood flow the river 
stream. readily seen from these definitions that the floodway must not 
obstructed its efficiency preserved. the other hand, the ef- 
ficiency the floodway may increased such measures enlarging the 
channel, removing obstructions, and other works. Such works should receive 
consideration some cases planning the use the flood plain and weighing 
the economic influences demanding the placing improvements the flood 
plain against the need restrict flood damages and preserve the degree 
efficiency the existing floodway. 

Such problems are highly technical and are usually too involved for the 
local zoning board solve satisfactorily. such situations, help should 
solicited from the local engineering department adequately staffed deal 
with the problem, from consultants, and from state federal water resources 
agencies. 

The action required individual cities and counties obtaining effective 
regulation flood plains has been very ably summarized 
points out that local communities need consider: 


Forming planning commission one does not already exist and pre- 
paring master plan land use the community. 

Reviewing existing ordinances and codes that contain provisions con- 
cerning development flood areas, revising them, and adding them they 
are inadequate regulate development under major flood conditions. 
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FLOOD PLAIN REGULATION 


Preparing zoning ordinances, subdivision regulations, and building 
codes, not already existence, that certain definite and adequate pro- 
visions for regulation flood-plain development are 

Cooperating with federal and state agencies the procurement 
needed flood data. 

Instituting program acquisition appropriate portions the most 
hazardous and frequently flooded lands. 

Instituting program annual maintenance stream channels. 


addition the actions cited above, the will something about the 
problem and the will enforce regulations are very necessary ingredients 
achieving effective results. 


State Legislation 


the discussion local regulation flood plains was pointed out that 
adequate statutory authorization for planning, zoning and subdivision control 
the city, county and regional levels lacking number states. Such 
authorization the state first step toward enabling local communities 
assume the duties that must carried out the local level. 

addition granting authority for local action, the states can take other 
steps support effective flood plain regulation. This support can the 
form promotion local responsibility flood plain regulation, technical 
and advisory help with flood plain problems, state supervision certain types 
developments, such dams and levees, and state control floodways 
under certain circumstances. 

Lack understanding flood plain problems, ignorance the potential 
flood danger the flood plain, and local pressures disregard recognized 
danger account for much the inaction dealing with flood plain problems. 
effective program contacting planning officials and zoning boards in- 
form them the flood dangers their communities, steps they may take 
combat those dangers, and encouraging them adopt the plans and or- 
dinances necessary can carried out the state level. state agency with 
one two people, expert these matters, can serve all the communites 
the state promoting better flood plain regulation. 

The determination the floodway requirements and boundaries tech- 
nical task that involves study past and probably future flood discharges, 
characteristics the floodway, and possible modifications the floodway. 
History past floods frequently insufficient furnish safe guide for 
floodway capacity requirements, and future greater floods must considered. 

The technical experience for making such studies will lacking most 
planning boards and most city engineering staffs. However, most states 
now have water resources agencies with responsibilities for making studies 
floods and other water problems. These agencies, because their famil- 
iarity with such problems can provide assistance and advice local officials. 

Local zoning regulations should prohibit encroachments floodways and 
these regulations should supported state encroachment laws. addition, 
the state encroachment laws should also protect those floodways not protected 
local laws. Local planning and zoning officials may find advantageous 
require the approval state agencies prior condition granting ex- 
ceptions floodway zoning ordinances. 
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State laws should prohibit encroachment floodways and should provide 
for adequate technical review and approval all structures, such bridges, 
dams, pipelines, etc., that must constructed floodways. All filling, ex- 
cavating, other constructions that might adversely affect the efficiency and 
capacity the floodway should prohibited legislation and penalties 
should provided for violations. 

Dams, which can have high hazard factor and can serious menace 
developments flood plains and public safety, lack adequate regulation 
many states. Only states require permits for the construction dams, 
require approval plans before construction, and only provide for inspec- 
tion after construction. Even though inspections are made, many states lack 
adequate statutory authority order the abatement removal hazardous 
conditions assess penalties for noncompliance with orders safety 
measures. 

Levee systems protecting urban areas frequently suffer from inadequate 
design and maintenance with consequent danger life and property. Also, 
such structures often lack responsible owner provide necessary mainten- 
ance. This situation can best remedied through state inspection and 
requirement that local communities provide adequate maintenance. When 
local maintenance deficient, the maintenance should performed the 
state and charged back the local community. 

Flood data can best gathered and analyzed the state level because 
full time specialized technical staffs, which cannot financed lower levels 
government, can provided for such purposes this level. addition, 
the restrictions local governmental boundaries are not limiting factors 
the extent and scope investigations. The controlling features that influence 
flood heights are frequently considerable distance beyond the urban 
county limits that limit local jurisdictions. 

would appear that legislation the state level should provide: 


Statutory authority for the creation local planning agencies with 
prepare comprehensive community plans including designation flood 
plain and floodway areas. 

Statutory authority for enacting zoning and subdivision control ordinan- 
ces. 

floodways. 


Enactment laws governing construction and maintenance dams and 
levees. 

Establishment state flood control water resources agency with 
regulatory powers, authority make investigations, collect basic data, pre- 
pare flood-evaluation reports, provide advisory service local planning 
agencies, and cooperate with federal agencies. 


The role the state should support the local governmental units 
any policing actions necessary for local regulation, and provide state regu- 
lation those matters that require uniformity action for justice and prac- 
tical operation and those situations that transcend the limits local 
ernmental jurisdictions. 
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FLOOD PLAIN REGULATION 


Federal Legislation 


The historic approach the federal government the problem reducing 
flood damages has been through the construction flood control flood pro- 
tection structures. Little effort has been exerted prevent the creation 
new flood damage areas encourage local responsibility this field. 
However, there awakening the federal level the need for preventing 
the creation new sources flood damage addition continuing the tra- 
ditional effort provide protection old areas. 

The Tennessee Valley Authority its recent report Program for Re- 
ducing the National Flood Damage Potential” cites the futility continuing 
provide structural measures for flood reduction flood protection certain 
areas and the same time permitting new damageable developments take 
place even faster rate other areas. 

The TVA report suggests three-pointed plan for placing more emphasis 
the control the growth the flood damage potential: 


expanded program for collection, analysis and wider use flood 
data. 

broad approach flood damage abatement which includes preven- 
tative well corrective measures. 


More active participation state and local governments solving 
the flood problem. 


implement the program the TVA recommends that: 


national flood damage prevention policy should developed. The 
objectives the policy should (a) prevent the unnecessary spread 
buildings and other improvements new flood hazard areas, combined with 
continued efforts alleviate damages those flood hazard areas which 
have already been developed, and (b) encourage appropriate state and 
local bodies assume the responsibility for initiating and developing plans 
and programs for local flood damage abatement through construction 
protective works and establishment land use controls necessary. The 
policy should also provide basis for technical and financial assistance 
federal agencies. 

protective works are economically justified part the local 
program for flood damage abatement, the federal financial contribution 
such projects should contingent upon the requirement that the city 
state initiate and prepare engineering plans the proposed works and that 
the plans submitted appropriate federal agency for review and 
approval. 

prerequisite for the contribution federal funds local flood 
protective works, state and local governments should required adopt 
and administer flood plain zoning and other such controls are needed 
prevent the unnecessary spread buildings and other improvements 
areas subject flood damage. 

Provision should made insure that all federal agencies having 
responsibility related site selection financing physical structures 
observe local flood plain zoning and other standards carrying out their 
responsibilities, and cities without land use controls but subject flood 
damage require establishment such controls before federal programs 
are approved and put into execution. 


expanded and systematic program collecting and analyzing 

basic flood information existing federal agencies should authorized 
which will meet the needs states, local communities, industries and 
others adjusting flood conditions. 

For the purpose providing cities and communities with basic flood 
information which prepare plans for growth and development and for 
the prevention local flood damages, appropriate federal agencies should 
authorized, where such authorization necessary, prepare upon re- 
quest local flood studies which would describe and analyze the local flood 
situation. 


The TVA recommendation, addition providing for greater federal ac- 
tivity collecting basic data and providing basic flood information local 
communities, would make flood plain regulation requirement for qualifying 
for certain types federal assistance. This requirement would provide con- 
siderable stimulus local regulation. 

The Chief Engineers, Corps Engineers, has recently instructed 
all the Division and District offices give greater consideration flood 
plain regulation. directed that survey reports the Corps Engineers 
will discuss adequately the problems encroachments channels and de- 
velopments flood plains. Where encroachments and developments can 
expected, appropriate practicable provisions are included proposed 
plans improvements and requirements local cooperation protect the 
flood-carrying capacity channels and floodways. 

These moves federal agencies are the right direction and general 
need legislation place them effect. However, there are deficiencies 
the general policy level that detract from their effectiveness. policy dec- 

laration Congress that would endorse the recommendations the TVA and 
the directive the Chief Engineers and extend the recommended practices 


all federal agencies would greatly assist promoting better flood plain 
regulation. 


SUMMARY 


Perhaps the greatest hindrance accomplishing more widespread flood 
plain regulation the lack adequate authority for planning and zoning 
city, county and regional levels many states. The acceptance planning 
and zoning methods for flood plain regulation still its infancy and 
greater educational effort needed bring about greater use these tools 
for reducing the flood damage potential. 

State laws provide adequate planning and zoning authority should en- 
acted promptly where necessary. Such laws should made adequate for the 
city, county and regional levels. 

Although all states have provisions for administering water resources ac- 
tivities some state agency some degree, many the agencies lack ade- 
quate powers deal with flood problems, collect basic data 
regulations. addition, many are inadequately financed and staffed. The 
states’ effort should suuport local flood plain regulations, collect basic 
flood data and prepare flood evaluation reports, provide advice and technical 
service local planning agencies flood plain matters, administer state 
and dam-safety provisions, provide levee inspection 
and require adequate levee maintenance for urban areas, and cooperate with 
federal agencies. 
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FLOOD PLAIN REGULATION 


The federal government should adopt policy that will require local regula- 


tion flood plains where indicated federal flood control projects, enlarge 
the program for the collection basic flood data, prepare flood-evaluation 

reports particularly regions that cross state boundaries, and establish cri- 
lood teria for determining the extent which loans private construction flood 
plain areas would guaranteed federal agencies. 
re- 
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IMPROVED TUNNEL SPILLWAY FLIP BUCKETS 


SYNOPSIS 


This paper discusses several buckets, the type used deflect flip 
tunnel spillway discharges downstream, terms their desirable unde- 
sirable features. Several new types buckets developed from hydraulic 
model tests are then described. Using data from these tests, dimensionless 
curves are presented aid determining the jet trajectory length, the 
spreading the jet, the tail water drawdown the bucket, and the pressures 
the floor and side walls the bucket. 


INTRODUCTION 


tunnel spillway composed two basic parts—an upstream spillway 
crest, free controlled, and downstream tunnel, part which sloping 
and part near horizontal. From the standpoint economy the tunnel diameter 
must kept minimum. Since the tunnel never allowed flow full be- 
cause the possibility siphonic action producing dangerous flow conditions, 
necessary keep flow velocities high and prevent turbulent areas 
the tunnel. Spillway tunnels are usually designed flow from 3/4 7/8 full 
maximum discharge, making the outflow the tunnel portal relatively deep. 
The combination depth and velocity produces the highest possible concen- 
tration energy and increases the difficulty obtaining satisfactory flow 
conditions where the flow spills into the river. example, the Glen 
Canyon tunnel spillways,! the maximum discharge 276,000 cfs produces 
159,000 horsepower per foot width the tunnel portals. Grand Coulee,2 


Note: Discussion open until May 1960. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. Paper 2316 
part the copyrighted Journal the Hydraulics Division, Proceedings the 
American Society Civil Engineers, Vol. 85, No. 12, December, 

Hydr. Engr., Div. Eng. Labs., Bureau Reclamation, Denver, Colo. 

Hydr. Engr., Head, Hydr. Structures Unit, Div. Eng. Labs., Bureau 
Reclamation, Denver, Colo. 
See Table for references 1-8. 


4 
q 


4 

q 


overfall spillway, where the maximum discharge 1,000,000 cfs, the ener- 
per foot width only 15,650 horsepower, one-tenth that Glen 
Canyon. 

were feasible construct efficient hydraulic jump stilling basin 
the end one the Glen Canyon tunnels, the basin depth, from apron tail 
water elevation, would need 170 feet deep. The hydraulic jump length 
would over 1,000 feet and would require basin 700 800 feet long 
more. Basin appurtenances such baffle piers could not used effectively 
because the high entrance velocity, 165 feet per second, would produce cavi- 
tation problems. The cost structure this size would prohibitive, and 
readily seen why other types structures are used the end tunnel 
spillways. Buckets have been the most common these structures and were 
probably derived from the slight upturns placed the base early overfall 
spillways. not clear whether the designers intended that these buckets 
operate free submerged. some cases, the upturn was too slight pro- 
duce measurable effect thick jet, but probably the intended purpose was 


deflect the jet downstream prevent undermining the spillway structures, 


Buckets this type are referred variously “ski-jump,” deflector, dif- 
fuser, trajectory, flip buckets. For ease expression, the term flip 
bucket will used this paper. 

Flip buckets are not substitute for energy dissipators because such 
bucket inherently incapable dissipating energy within itself. The purpose 
flip bucket throw the water downstream where the river-bed damage, 
which usually certain occur, does not endanger the safety the dam, 
powerplant, other structures including the flip bucket itself. accomplish- 
ing this primary function, buckets are also designed spread the flow across 
much the downstream channel considered desirable order re- 
duce riverbed damage much possible. The jet trajectory modified 
necessary cause the jet impinge the tail water surface the desired 
location, and when possible, the steepness the jet trajectory the point 
impingement selected produce horizontal and vertical velocity compo- 
nents which produce most favorable flow conditions the river channel. 

Although with the present state knowledge impractical generalize 
the design flip buckets, intended this paper present certain basic 
facts which have been found true result extensive hydraulic model 


Table 1 
DESCRIPTION OF PROJECTS 
Reference: Name and : Maximum headwater : Tunnel 
No. location Agency : @ischarge : to dimensions 
bucket invert: 


1 :Glen Canyon Dam :Bureau of Reclamation : 276,000 cfs: 588" :2 tunnels 
: Colorado River Storage Hy : 3 : each 41' diameter 


:Grand Coulee 2 pillway 
: Columbia Basin Project : 1,650" wide 


3 :Hungry Horse ;Bureau of Reclamation : 50,000 cfs: =) 731° diameter 
: Hungry Horse Dam Project 
Montana 


: Miseouri River Basin Project: : 3 : horseshoe conduit 


: :2 horseshoe tunnels 
: Ministry of Agriculture : : : 37.08' diameter 


:Taiwan Power Company 3 :27° diameter 


: North Carolina : : : : each 34' diameter 
:Trinity Dem 2 00 5 diameter 


: Central Valley Project 
: California 


of 
> ~ 
B 
; Project, Arizona 3 C2 
TYellowtail Dam TBureau of Reclamation 173,000 cfs: Sie" diameter 
7 :Fontana Dam Tennessee Valley Authority 3 , 000 cis: tunnels 
4 


testing and prototype observation. hoped that extension the ideas pre- 
sented, the form discussions, will help provide better understanding 
the requirements necessary for improved flip bucket design. 


Bucket Design Problems 


usually difficult impossible predict the flow pattern expect- 
from particular bucket mere inspection the bucket shape. Because 
variations velocity and depth, the spreading and trajectory character- 
istics given, bucket can determined only testing hydraulic model. 
Since the authors have had the opportunity test various types buckets and 
observe hear first hand their performance the field, the findings 
these tests should, therefore, interest designers who must often se- 
lect bucket type before the hydraulic model tests are made. 

the course developing and improving bucket designs, number 
difficulties have been found and overcome. The following examples indicate 
the problems which may encountered bucket design and which may not 
generally known. 

The flip buckets the tunnel spillways Hungry Horse and Yellow- 
tail the Bureau Reclamation projects, and Yanhee and Wu- 
Sheh Dam” being built Thailand and Formosa, respectively, are similar and 
are what may called “standard” type. The buckets are placed down- 
stream from transition which changes the circular horseshoe-shaped 
tunnel flat bottom correspond the flat bottom the bucket. High ve- 
locity flow the tunnel makes difficult design short transition; long 
transitions are usually costly. the transition not carefully designed, and 
preferably checked model studies, there the possibility dangerous sub- 
atmospheric pressures occurring the corners. The transition, therefore, 
becomes much design problem the bucket. 

The Fontana spillway buckets not have upstream transition. 
The bucket inverts are circular, the same the tunnel inverts, Fig. 1.* The 
buckets were shaped trial 1:100-scale model tested the TVA Hy- 
draulic Laboratory. The curved surfaces the finally developed buckets 
could not defined ordinary dimensioning even mathematical 
equations. That the buckets were well designed has been proved subsequent 
operation the structure, but the methods necessary convert the model di- 
mensions scale 1:100 prototype dimensions were quite laborious. 
Because the high velocity flow the bucket, dimensions taken from the 
model could not “scaled up” directly. Any small irregularity misaline- 
ment when multiplied 100 could have been sufficiently large produce 
cavitation the prototype bucket. was, therefore, necessary convert the 
dimensions 1:10-scale bucket, and after smoothing these, convert the 
corrected dimensions 1:1 scale. 

some buckets, particularly those foreign dams, serrated toothed 
edge has been placed the downstream end the bucket. The teeth are 
provide greater dispersion the jet before strikes the tail water surface. 
High velocity flow passing over the sharp edges may produce cavitation 
damage the concrete surfaces. 


*Photographs from Technical Monograph No. 68, Hydraulic model studies, 
Fontana Project, Tennessee Valley Authority. 
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The problem draining tunnel which has flip bucket the downstream 

end provides challenge design. The drain must placed surface ex- 
posed high velocity flow. Even though possible design develop 
the laboratory drain opening which will not produce cavitation pressures, 
difficult obtain field construction the necessary tolerances prevent 
cavitation from occurring. ideal bucket design would self-draining and 
would not present cavitation problem the drain structure. 


Improved Bucket Designs 


Recently number tunnel spillway flip buckets have been developed 
the Bureau Reclamation Hydraulic Laboratory that seem offer simple but 
very effective methods directing the flow away from the structure and which 
also overcome, part, the difficulties described the preceding paragraphs. 
Although single bucket eliminates all the undesirable features 
bucket, the use the principles described the following pages will help the 
designer provide improved bucket particular structure. review, 
ideal bucket should provide (1) easy drainage the tunnel, (2) bucket 
shape which can defined and expressed prototype size ordinary di- 
mensioning ordinary drawings, (3) need for upstream transition, and 
(4) impingement area which may shaped, simple additions basic 
bucket, fit the existing topographic conditions. Some the buckets de- 
scribed are unique and probably cannot generally used without some adapt- 
ation; however, the others are basic type and have only minor additions 
accomplish some specific function. 

One the unique designs was the Trinity Dam8 spillway bucket developed 
1:80-scale model. The spillway tunnel enters one side wide shallow 
river channel and the flow tends cross the river diagonally. was neces- 
sary discharge the flow into this channel without creating excessive eddies 
that might erode the riverbanks cause disturbances the vicinity the 
powerhouse tailrace. The spillway uncontrolled morning-glory and, 
consequently, the flow can vary from few second-feet maximum 
22,000 cfs. The velocity the bucket 122 feet per second. Because small 
flows may occur for days, was desirable for low flows leave the bucket 
close the riverbed elevation possible prevent excessive erosion 
near the base the structure. the other hand, large flows should 
flipped downstream away from the with much dispersion possi- 
ble prevent erosion and induced eddies from damaging the structure. the 
usual flip bucket, hydraulic jump forms the bucket for small flows and the 
water dribbles over the bucket end and falls onto the riverbed. This could 
cause erosion which would undermine the structure. When the jump first 
swept out the bucket, the jet usually lands near the structure and erosion 
and undermining the structure may still occur. Trinity Dam, the foun- 
dation conditions the end the tunnel were such that was deemed neces- 
protect against the possibility erosion and undermining. order 
place the bucket near riverbed level, the semicircular channel constructed 
downstream from the tunnel portal was curved downward trajectory curve, 
and the flip bucket structure was placed the end, Fig. The flip bucket 
surface consisted three plane surfaces placed that they spread and shape 
the jet fit the surrounding topography. Large flows are spread into thin 
sheet having contact line with the tail water surface considerable distance 
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downstream, Fig. 3A. However, even small flows are thrown downstream 
well away from the base the bucket. 

training wall was used prevent spreading the jet the high, land 
side, the bucket. There was wall the low river side the bucket. 
flows less than 1,000 cfs, hydraulic jump formed over the horizontal 
surface and part way the slope the bucket; the flow spilled out the low 
side the bucket into the river channel. The open side the bucket was only 
feet above the river. Had the flow been confined both sides and 
forced spill out the end, the drop would have been over feet and addition- 
protection the bucket foundation would have been required. discharges 
greater than 1,000 cfs, the jump swept out the bucket without hesitation and 
with sufficient velocity that the flow was carried well downstream away from 
the structure. the discharge increased, the jet was flipped farther down- 
stream and became increasingly dispersed. The long contact line between the 
jet and the tail water reduced the unit forces the tail water, and the eddies 
induced the ends the contact line were thereby found minimum. 
Since one side the bucket entirely open, the bucket self-draining. 
Other advantages this design are that the bucket may defined for proto- 
type construction with few simple dimensions, and curved warped 
forms are necessary for prototype construction. 

Another unusual type flip bucket was developed for the Wu-Sheh Dam 
tunnel spillway. Construction schedules and geologic conditions the field 
made necessary modify this bucket from the standard type described 
earlier this paper. After the line the tunnel had been established and 
construction the tunnel started, was found necessary result model 
tests change the direction the flow entering the river channel. Earth and 
rock slides, during the diversion period, made necessary construct re- 
taining walls the tunnel portal area which restricted the length the flip 
bucket. Hydraulic model studies were made determine how much turning 
the jet was required and whether the turning could accomplished the 
tunnel. The tests showed that was undesirable turn the tunnel and that all 


A> 
{ 
\ 
\ ~ 
\\ 
\ \ \\ 


DISPERSION FLIP BUCKET 


4 
wee 
: 
= 
~ 
Flow 
S& 
SS SS 
\ 
| = 
| 
aa 
\\\\ 
Fig. 


FLIP BUCKETS 


13,000 cfs, Basin Acts Flip Bucket 


Fig. 3B. Combination Hydraulic Jump Basin-Flip Bucket 
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turning should accomplished the bucket. The final bucket determined 
from model studies used curved turn the flow—a batter the left 
wall prevent congestion the bucket and reduce hydraulic loads the 
larger discharges, and fillet the junction the left wall and floor 
smooth control the jet undernappe, Fig. The resulting bucket was 

“tailor-made” direct the flow impinge near the middle the river 
channel and obtain the greatest dispersion possible all discharges. The 

surfaces this bucket could also defined ordinary dimensioning. 

Piezometers placed the side walls the bucket showed 
atmospheric pressures all discharges. The maximum pressure recorded 
the left wall was feet water, Fig. Before the wall was battered, the 
maximum pressure probably would have been much larger due more di- 
rect impact the converging wall. 

The Yellowtail Dam tunnel spillway flip bucket dual purpose bucket 
similar some respects the standard buckets. The tunnel curved 
bottom horseshoe-type conduit. Two hundred fifty feet upstream from the 
portal, the tunnel changes flat bottom horseshoe conduit, and the invert 
drops feet means combination transition-trajectory curve 172-1/2 
feet long. The bucket has flat horizontal floor 130 feet long and 62-1/2- 
foot-long upward sloping sill, Fig. spillway flows 12,000 cfs, 
hydraulic jump forms the bucket and relatively quiet water discharged 
into the downstream channel. the spillway discharge increases, the jump 

moves downstream and 13,000 cfs sweeps out the basin, Fig. 3B; for 
greater discharges and the maximum, 173,000 cfs, the basin acts 
flip bucket. The basin bucket placed low solid rock that discharges 
the unstable zone, 12,000 13,000 cfs, cannot undermine the structure. 
This basin was developed the Hydraulic Laboratory serve the specific 
purpose acting hydraulic jump basin for the most prevalent spillway 
discharges—discharges expected exceeded only every years, and act- 
ing flip bucket prevent damage the structures during large floods. 
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The reason for using the hydraulic jump for part the discharge range was 
protect the river channel against clogging with talus which was present 
the canyon large quantities and was expected move high velocity 
stream contacted it. Following large discharges, was expected that re- 

opening the channel achieve full power head would necessary. 

The flip buckets for the Glen Canyon Dam tunnel spillways are example 
buckets developed eliminate the tunnel transition and the need for flat 
bucket invert. The buckets the portals the 41-foot-diameter tunnels are 
opposite sides the river and are aimed discharge acute angles with 
the center the river. The left bucket farther downstream than the right. 
Each bucket designed handle the maximum discharge 138,000 cfs ata 
velocity about 165 feet per second. This represents over 13,000,000 horse- 
power energy released into the river during maximum discharge. 

the preliminary design, there was 70-foot-long transition between the 
circular tunnel and the rectangular channel containing the flip bucket. Hy- 
draulic model studies indicated that the transition was too short, and that sub- 
atmospheric pressures would sufficiently low produce cavitation and 
damage the structure. Two alternatives were developed during the model 
studies; one was use 100-foot-long transition which the change cross 
section was accomplished without dangerous pressures occurring and the 
other was eliminate the transition continuing the circular tunnel invert 
downstream intersect the upward curve the flip bucket. The latter 
scheme was developed and will used the prototype structure; identical 
buckets will constructed the twin spillways. effect, the transition and 
the bucket are combined into the bucket structure without complicating the de- 


sign the bucket. 

Because the flat-bottom portion these buckets diverges plan, small 
flows are spread laterally more than for the flat-bottomed bucket. the dis- 
charge increases, the rate spreading decreases that easier ac- 
commodate the jet for flood flows relatively narrow channel. Fig. 
shows comparison the flow from the two types buckets. the flat- 
bottomed bucket which preceded transition, the flip curve extends 
across the full width the bucket for its entire length; all the flow 
given elevation are turned simultaneously. the alternate bucket, the 
flip curve turns the lower flow elements the center the stream first and 
gradually widens its zone influence the flow moves downstream, result- 
ing greater dispersion the jet. effect, the flow along the center line 
the bucket turned upward while the flow elements either side the 
center are turned upward and laterally. Training walls may used limit 
the lateral spreading. subsequent testing, deflectors were added the 
training walls make the jets conform the shape the river 
channel and surrounding topography, Fig. 

The flip bucket used the Flaming Gorge Dam tunnel spillway was the 
type used the Glen Canyon spillways. The maximum design flow 
for Flaming Gorge spillway 28,800 cfs; the velocity the flow the portal 

the 18-foot-diameter tunnel about 140 feet per second. The energy the 

jet the flip bucket equivalent 1,000,000 horsepower. operation the 

flow appearance the Flaming Gorge bucket was entirely different than that 

the Glen Canyon buckets. The Flaming Gorge jet was well dispersed the 

lower discharges and became more compact the discharge increased, Fig. 

The Glen Canyon jets were well dispersed for all flows, and the change 
lateral spreading with discharge was not apparent. the Flaming Gorge 
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Transition Flip Bucket with Side Wall Deflector 
Flow Froude Number 5.64 


Typical Jet Profile 
35° Transition Flip Bucket 


Fig. Glen Canyon Dam 
Flip Bucket Studies 
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FLIP BUCKETS 


Jet Well Dispersed Low Flows. 
Froude Number 10.3 


More Compact Jet Forms for Larger Flows. 
Froude Number 6.8—Maximum Discharge 


Fig. Flaming Gorge Dam 
Flip Bucket Studies 
35° Transition Bucket 


bucket, the water rose the sides the bucket low flows, forming ef- 
fect “U”-shaped sheet water which the bottom and sides were equal 
thickness. The vertical sides the “U” followed the line the bucket side 
walls after leaving the bucket, while the bottom sheet water had tendency 
diverge either side. The vertical fins had shorter trajectory than the 


lower sheet and falling would penetrate the lower jet, tending spread 
disperse it; this can seen the photographs Fig. the discharge 
increased, the size the fins relative the thickness the lower sheet be- 
came insignificant and longer had this spreading effect. The differences 
the Glen Canyon and Flaming Gorge jets might explained the fact that 
the flow depth for maximum discharge was about per cent the diameter 
the Glen Canyon tunnel and per cent the diameter the Flaming che 
Gorge tunnel. For flow 0.61D Flaming Gorge, the jet was still well dis- flo 
persed. 
Both the Flaming Gorge and Glen Canyon buckets were modified reduc- 
ing the height the river side wall. The Flaming Gorge bucket located ret 
well above the maximum tail water elevation that the wall could cut down 
the spring line the tunnel invert curve without tail water interference. buc 
The effect was eliminate the fin that formerly rose along the wall. The jet use 
spread out evenly the right and was better dispersed than before. The Glen 
Canyon buckets are located more closely the maximum tail water elevation nur 
and prevent the tail water from interfering with the jet, the river wall 
cut down only feet above the spring line the tunnel invert. Sufficient 
wall remained train the jet and very little difference the flow pattern ver 
could detected. the 
div 
Design Considerations 
sig 
Tunnel spillways usually make use part the river diversion tunnel. con 
The downstream portion the diversion tunnel becomes the horizontal portion 
the spillway tunnel—the bucket added after diversion needs have been 
satisfied. Since the diversion tunnel one the first items construction, 
often necessary, because time limitations and construction schedules, 
determine line and grade for the diversion tunnels before the details the 
spillway are known. Care selecting the exact position and elevation the and 
diversion tunnel, while keeping mind its ultimate use spillway tunnel, 
will help provide dual purpose tunnel which will satisfy the temporary sho 
well the final demands with the least amount modification when the sar 
bucket added. tra 
The following sections cover the items which should considered during 
design and which will help provide simple bucket structure having desira- the 
ble performance characteristics. 
tes 
Elevation Bucket Invert 
desirable construct the bucket and tunnel inverts the same ele- abl 
vation. Since diversion requirements make necessary keep the wid 
tunnel low provide the diversion capacity, the greatest danger that the 
tunnel will set too low for ideal spillway operation. This will require build} 
ing the bucket lip prevent the tail water from submerging the bucket. lim 
general rule, maximum tail water should higher than the elevation 


the center line the tunnel. the bucket set lower, difficulty may 
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experienced obtaining free flow low spillway discharges. The shape 
the tail water curve will determine the exact requirements. The drawdown 
tail water elevation the bucket caused the ejector action the jet may 
also affect the vertical placement. Drawdown discussed later. 


Flow Direction 


The bucket center line should continuation the tunnel center line, 
and the portion the diversion tunnel used for the spillway tunnel should 
straight. is, therefore, desirable aim the diversion tunnel that may 
used without change for the spillway tunnel. The tunnel direction should 
set that spillway flows will aimed downriver and that the design dis- 
charge impinges the tail water the center the discharge channel. The 
flow should directed minimize the diameter induced eddies the sides 
the jet since these can very damaging channel banks. ideal ar- 
rangement, the jet will wide the channel that there will little 
return flow from the downstream tail water. 

Fig. shows the angle divergence one side the jet leaving the 
bucket for two types buckets—the flat-bottom type and the transition bucket 
used Glen Canyon and Flaming Gorge spillways. both cases, the angle 
divergence plotted versus the angle inclination for range Froude 
numbers (of the flow entering the bucket). The flat bottom bucket produced 
very little change angle divergence for range Froude numbers 
inclination angles. The transition bucket showed considerable change di- 
vergence angle—from 12° for Froude number range 11. Since 
the higher Froude numbers occur low discharges, the transition-bucket-jet 
divergence greatest low flows. the discharge increases, the Froude 
number becomes smaller and the divergence angle decreases. most de- 
signs this favorable characteristic and results improved river flow 
conditions for all discharges. 


Drawdown 


For the conditions just described, the jet will act ejector lower the 
tail water upstream from the jet impingement area. From the Hungry Horse 
Dam model tests, feet drawdown was predicted for 35,000 cfs discharge, 
and was recommended that weir constructed the powerplant tailrace 
prevent unwatering the turbines. Prototype tests made for 30,000 cfs 
showed feet drawdown and demonstrated that the weir was indeed neces- 
sary. Hungry Horse the flow leaves the bucket 15° angle, making the 
trajectory relatively flat, Fig. 10; the jet wide the downstream 
channel. The drawdown maximum under these conditions. Glen Canyon 
the spillway jets not occupy the entire width channel, but the jet trajecto- 
steeper, and the discharge considerably greater. Hydraulic model 
tests have indicated that feet drawdown may expected. 

Other hydraulic model bucket tests have shown the drawdown appreci- 
able, particularly when the jet occupies large proportion the channel 
width. means have been found calculate the amount drawdown 


expected, except making careful measurements hydraulic model. How- 


ever, using measurements obtained several model studies and from 
limited prototype observations, the curve Fig. was derived. pre- 
sented here means estimating the drawdown that can expected with 
tunnel spillway and flip bucket. 
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Hungry Horse Spillway Tests. Spillway Discharge 
30,000 cfs. Side View Jet. 


Hungry Horse Spillway Model—Discharge 35,000 cfs 
Model-prototype Comparison 


Fig. 10. Hungry Horse Spillway 
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The intensity the ejector action and the resultant lowering the tail 
water the bucket have been found function the energy the jet 
and the amount resistance encountered when the jet strikes the tail water. 
the curve Fig. the abscissa the cross-sectional area the river 
flow near the point impact the jet divided the cross-sectional area 
the flow the tunnel portal. The river flow area the product the differ- 
‘ence between the no-flow tail water elevation and the tail water elevation for 
the discharge being investigated, and the average width the river near the 
point impact. The area the flow the portal obtained dividing the 
spillway discharge quantity the average velocity. The ordinate the ratio} 
the amount drawdown the depth tail water. The depth tail water 
the same depth used determine the river cross-sectional area. 

The curve, defined the test points shown, indicates with reasonable accu- 
racy the drawdown each dam site for which data were available. The test 
points include various shapes and depths channel and various types 
bucket jets. Further, the two prototype tests Fontana and Hungry 
Dams showed good agreement between model and prototype test results. 
ever, predicting drawdown future sites, the curve should used with 
caution until more data are available. 

tu 
Effect Trajectory Shape 

addition the effects drawdown explained above, the jet trajectory 
important other ways. The angle the bucket lip with respect horizontal 
determines the distance the water will thrown downstream. However, the 
steeper the angle, the more the jet will broken and slowed down air 
resistance. Both these effects cause the jet enter the tail water 
steeper angle. With steep entry, the vertical component velocity will 
greater, and the jet will tend dig into the channel bottom. With flatter tra- 
jectories, the horizontal component will greater, the drawdown will 
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greater, and the forward velocity will higher. High velocity channel flow 
may persist downstream from the impingement area for considerable 
distance the channel bottom does not erode produce deep pool. High ve- 
locity flow along the channel banks will then occur. the bottom erodes, 
energy dissipating pool will formed, and flow downstream will smoother. 
Bucket flip angles are usually constructed from between 15° and 35°. Angles 
less than 15° not give enough lift clear the bucket structure, and little 
usually gained, from any viewpoint, increasing the angle beyond 35°. 

Fig. contains family curves which may used estimate the tra- 
jectory length for inclination angles 45° and velocities 160 feet per 
second. curves were obtained from the simple equations for the path 
projectile, sin For given angle the equation may sim- 
plified shown the equations the right the trajectory curves. For 
bucket entrance. estimate the curve the lower right Fig. may 
used. Here, H', expressed percentage the total head plotted 
versus the percentage maximum tunnel discharge. seen vary from 
about per cent for per cent maximum discharge, about per cent 
for maximum discharge. Maximum discharge considered occur when the 
tunnel about three-fourths full the outlet portal. The points which de- 
termine the curve have ratios vertical drop horizontal tunnel length, 
H/L, from 0.15 1.9. 

Trajectory lengths taken from these curves have been found reason- 
ably accurate when checked hydraulic models. Some difference between 
model and prototype trajectory lengths may expected occur, however. 
Little known regarding model and prototype trajectory length agreement, but 
from measurements estimated scaled from photographs, and from actual 
measurements reported the author, (1) appears that the differences are 
usually not critical nature. The prototype trajectory shorter than the 
model theoretical jet and has steeper angle entry into the tail water. 
The difference believed caused the greater air resistance en- 
countered the high velocity prototype jet. From sketchy information 
few structures, the trajectory length the prototype for per cent maxi- 
mum discharge believed per cent shorter than the model, 

Fig. 12. There also are indications that the difference becomes less the 
prototype discharge increases. 

determining the radius the bucket curve, necessary provide 
radius least four times great the maximum depth flow. This pro- 
vides incline sufficiently long turn most the water before leaves the 
bucket and provides assurance that the jet will thrown into the desired area 
downstream. 


Pressures the Transition Bucket 


Because the simplicity and effectiveness the transition bucket, will 
probably used many future tunnel spillways. Extensive pressure 
measurements were therefore made several buckets having two different 
inclination angles, 15° and 35°, indicate that the buckets were safe against 
cavitation pressures and provide data for structural design. The results 
these tests have been summarized Figs. 13, 14, and and may 
helpful making preliminary designs. 
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FLIP BUCKETS 
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F=68 


Measured pressure 
Theoretical pressure; 62.5) 
where 
Developed distance from PC.to piezometer 
Developed distance from end bucket 


Froude number, computed from and PC. 
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Fig. shows pressures along the center line the transition bucket 
floor. The envelope curve includes inclination angles from 15° 35° and 
flows the Froude number range 6.8 10.3, the usual range operating 


conditions. The maximum pressure was found slightly greater than giv- 


Gumensky from theoretical considerations. The theoretical pressure 
Pt, expressed: 


where 


This maximum pressure occurred about 0.6 the bucket length from the up- 


stream end. Pressures rapidly became less toward the downstream end 
the bucket and reached atmospheric the bucket lip. 


w= 


ANGLE 


Measured pressure end bucket 
Theoretical pressure (See figure 


bucket 

sidewall 
Bucket 


PRESSURES END BUCKET 
Fig. 
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For some tests piezometer placed just upstream from the bucket lip, 
Fig. 14, indicated below atmospheric pressures, phenomenon which has not 
been satisfactorily explained. Experiments model buckets showed that the 
pressure this piezometer was affected the shape angle the down- 
stream portion the bucket lip. The curve Fig. shows the relation be- 
tween pressure and the angle The curve indicates that for given angle 
inclination should 35° more insure atmospheric pressures 
above the lip piezometer. The curve also indicates that the pres- 
sure will atmospheric. This not practical solution, however, since 
the piezometer will then upstream from the lip and new problem 
will created the end the extended bucket. should noted that the 
bucket side walls extend beyond the lip piezometer shown Fig. 14. 

The curves Fig. indicate the pressures expected the side 
walls the transition bucket from the base the wall the water surface. 
For inclination angle 35° the maximum pressure about times 
great hydrostatic and occurs near the base the wall about the three- 
quarters point, x/1 0.75, the bucket length. the end the bucket, 

x/1 0.99 the maximum pressure only four times great hydrostatic. 
For 15° the maximum pressure four times greater than hydrostatic 
x/1 0.26, 0.55 and 0.80 and only twice great static x/1 0.99. 
Other data are presented for different bucket radii, R/1 values, and stations 
along the bucket, x/1 values. Although the data are not complete, sufficient 
information presented make preliminary structural design. the 
Flaming Gorge Dam spillway bucket, one side wall was cut down the spring 
line the tunnel without objectionable spreading the jet occurring when the 
flow depth exceeded the height the wall. This procedure simplified the 


structural design the bucket reducing the overall load wall which 
had rock behind it. 
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TIDAL CHARACTERISTICS FROM HARMONIC CONSTANTS 


Bernard 


ABSTRACT 


| 

Harmonic analysis tide observations determines harmonic constants 

(amplitudes and epochs) the harmonic constituents the tide. Formulas 
for using the constants obtain various ranges are furnished and techniques 
for obtaining typical curves are illustrated for various ports. These are 
compared with predicted curves. 


INTRODUC TION 


Many engineers have interest tides since tide observa- 
tions provide basic control hydrography, are important relative main- 
tenance rivers and harbors, and are factor pollution studies. re- 
cent years, projects involving storm-water levels and tidal-power develop- 

ment have emphasized the importance tides engineers. The tide tables 

furnish adequate predictions for most engineering problems but there are 
some aspects tides that can best generalized directly from the harmonic 
constants. 

Harmonic constants are the amplitudes and epochs the harmonic con- 
stituents that make the tide any place. Through modified Fourier 
analysis, using the periodicies the movements the moon and sun and sea- 
sonal meteorological variations, the observed hourly heights the tide for 
period observations are operated separate the tide into elementary 

harmonic constituents. 

The determination mean ranges, spring ranges (average near the times 
new and full moon), perigean ranges (average when the moon closest 


Note: Discussion open until May 1960. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. Paper 2317 

part the copyrighted Journal the Hydraulics Division, Proceedings the 

American Society Civil Engineers, Vol. 85, No. 12, December, 1959. 


Chf. Currents and Oceanography Branch, Tides and Currents Div., Coast 
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the earth), etc., and the typical shape tide curve are best determined from 
the harmonic constants. the purpose this paper indicate what har- 
monic constants represent and ways which they can used. 


General Description 


The tide the periodic rise and fall the water that results from the 
gravitational attraction the moon and sun acting upon the rotating earth. 
The tide-producing force, varying with the constantly changing positions 
both sun and moon relative the earth, can resolved into number 
simple periodic forces. simulate the complicated motions the moon and 
sun, there are substituted number hypothetical tide-producing bodies, 
each having fixed period for its circular orbit around the earth the plane 
the equator. Each hypothetical body produces constituent tide with 
period that can calculated accurately from astronomical data. Therefore 
the actual tide any place may conceived being made 
constituent tides, the period speed each being determined its ap- 
propriate hypothetical tide-producing body. 

Fig. shows the theoretical tide earth completely covered deep 
layer water when the moon the equator. Friction and inertia have 
been disregarded assumed that the waters have responded instantly 
the attraction the moon. Note that there high water bulge the side 
the attracting body, another high water bulge the opposite side, and 
low water belt circling the earth between, going from pole pole. The ex- 
istence bulge the opposite side can shown mathematically, 
the premise that the tide producing force any point the difference 
the gravitational attraction the tide-producing body particle that 
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HARMONIC CONSTANTS 


point and the attraction particle the center the earth. this mathe- 
matical development also shown that the tide-producing force varies in- 
versely the cube the distance, rather than the square the distance 
the gravitational attraction. 


Period Tidal Components 


Inasmuch the average interval between consecutive transits the moon 


(upper lower) hours and minutes, this the period the princi- 


pal lunar constituent. The component designated the symbol M2, the 


showing related the moon’s motion and the subscript showing that 
semidiurnal (period about half day). The principal solar constituent 
component named with period exactly hours. 

The components and will phase both full moon and new moon. 
From Fig. apparent that will make little difference whether the moon 
and sun are line the same side the earth (new moon) line with 
the center the earth but opposite sides (full moon), The large ranges 
new and full moon are called spring tides. With the moon first last 
quarter, the high water coincides with the low water. This results 
smaller than average semidaily ranges tide called neap tides. 

The periods and 12.42 and hours respectively (corresponding 


hourly angular speeds 28.98° and 30.00° per hour) will bring these com- 


ponents conjunction opposition the proper phases the moon. 

shown Fig. for Havre, France. Note that spring tides, with large 
ranges, not only are high waters higher than average but low waters are lower 


than average. neap tides, high waters are lower than average and low wa- 
ters higher. 


TIDE CURVE 
FOR HAVRE, FRANCE 


FIGURE 


The moon moves elliptical orbit rather than circular one. Since its ang 
distance from the earth varies, its tide producing force varies accordingly, per 
being greatest when the moon closest the earth (221,000 miles the 
and smallest when maximum distance from the earth (253,000 the 
apogee). The period from perigee next perigee averages about 27-1/2 
days. simulate this changing attraction the moon, introduce tha 
ponent with period 12.66 hours (speed 28.44° per hour) that will 
known the larger lunar elliptic semidiurnal component. Det 
The distance the earth from the sun also varies during the year although 
the degree variation considerably less than the moon from earth. Com- dee 
ponents representing this variation well numerous other lesser varia- Hov 
tions the apparent paths the moon and sun are used tide predictions 
but will not defined this paper. 
Thus far have dealt with semidiurnal tidal variations, restricting the the 
moon and sun the plane the equator. However the declination the sun tud 
varies annually between about 23.452°S and 23.452°N while the moon has ex- 
treme declinations 28.597°S and 28.597°N every years. 
Note Fig. what happens when the moon extreme north declination. 
The high water bulges are still present under the moon and the opposite 
side. The low water belt still circles the earth between, but longer 
pole pole. When the earth rotates its axis, point longer experiences 
the same high tide with both upper and lower transits the moon. some 
latitudes there are two high waters unequal height and nother latitudes 
just one high and one low water day. Thus there diurnal variation de- 
pending the declination the tide producing body. The principal diurnal 


combined luni-solar diurnal component with period 23.93 hours and 
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angular speed 15.04° per hour. lunar diurnal component with 
period 25.82 hours and angular speed 13.94° per hour. The sum 
the and angular speeds equals the angular speed, that each day 


the diurnal high water comes the same position relative the two semi- 


diurnal high waters. Furthermore the speeds and are chosen 
that they are conjunction when the moon extreme north south dec- 
lination and opposition when the moon the equator. 
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Determination Tidal Characteristics 


the early part this paper there was specified earth with uniformly 
deep water and resulting instantaneous response tide-producing body. 
However our earth does have land masses impeding the progression tide 


waves and the water areas are not uniformly deep. consequence the 


theoretical development gives only the basic periods the components 
the tide and same approximate relations between components but the ampli- 

tudes and epochs must determined from observations special form 
harmonic analysis. Fig. shows complete period component tide. 

The amplitude the average semi-range, CF, the vertical height from 


tical tide-producing body, the time from maximum high water the 
Three systems reference for the epoch component are used, each 


which designed for particular purpose. Since the available constants 
may referred any these, each will briefly mentioned. repre- 


the time passage the hypothetical tide-producing body over the 
the place for which tidal constants are furnished, the epoch 
marked for Kappa. represents the time passage over the 
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FIGURE 
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meridian Greenwich, the epoch noted for Greenwich transit. Fi- 
nally represents the time passage over the meridian the local standJ 
ard time zone (for example 75° for Eastern Standard time New York), 
the epoch noted lower case “g” and the corresponding time zone must 
indicated. Formulas for computing either given any one the the 
three, are furnished Coast and Geodetic Survey publications. 

Once amplitude and epoch are determined for each component, each one} 
can represented cosine curve. These curves are summed 
tide-predicting machine and tide tables can prepared for many years 


advance provided that the related hydrography water basins does 
change. Ordinarily about components are used for station the Coast con 
and Geodetic Survey although the machine designed for components. 


Various combinations components are used for different stations and all ious 
components significant amplitude, obtained from analysis tide Coa 
vations, are used for each station predicted. star 
Examples Combining Various Components 


The amplitudes and epochs the five components described previously, 
M2, S2, N2, Kj, and can combined approximate typical shape 
curve, and approximate various ranges and mean intervals. 

noted previously the combined hourly speed and 
hourly speed Consequently, although the combined and curve 
will vary amplitude depending the declination the moon, their combin 
high water will always fixed position relative the two semidaily high 
water Mg. Fig. shows some possible combinations and 
curves given equal amplitudes and various phase relationships the epochs, 
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FIGURE 
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the same meridian. Note that when the inequality 

the high waters, the low waters being equal. When the difference 

rk), 90°, the inequalities are equal and the sequence higher high lower low. 
must with difference 180°, the high waters are equal, the inequality being 

the the low waters. With difference 270°, the inequalities are again equal, 

but the sequence lower low higher high. Using these four sets curves 
ich guides, possible visualize the relationship for any value Mg° 
combining the epochs, multiples 360° may added arrive 
years value between and 360°. Similar relationships will hold for and 


not amplitudes that are unequal; Fig. may aid visualizing various 
combinations. 

ts. amplitudes the harmonic constants can used approximate var- 
all ious average ranges. These have been derived from formulas and tables 


Coast and Geodetic Survey Special Publication 260, Manual Harmonic Con- 
Reductions, 1952. 


Mean range (long period average) 

Spring range (at new and full moon) 2.1 
Neap range (at quadrature) 2.1 (Mg 
peof Perigean range (moon closest earth) 2.2 

Apogean range (moon farthest from earth) 
Tropic range (diurnal range extreme decli- 
ind the epochs the various components that combine form spring tides, 
tides, etc., were equal, these particular ranges would occur exactly 


new and full moon, perigee, etc. This not ordinarily the case, and the 

lag can calculated exactly from the epochs. The lag hours for spring 

jand neap tides found the formular 0.984 the lag perigean 
apogean tides, the formula 1.837 and the lag tropic 
0.911 positive lag means that the tide follows the ap- 

astronomical event; negative lag indicates that the tide precedes 
event. 


Some tidal components significance are not astronomical nature. 
are harmonics Mo, their angular speeds being exactly two and three 
that These are called shallow water components and reflect the 
the cosine-type tide curve because the physical features 
basins. The period basin determined its length and depth. 
the period basin approximates harmonic the tidal period, in- 
repetitive disturbing wave that more less uniformly distorts the 
the tidal wave. Figs. and show the shape the curve affected 
and respectively. The degree distortion depends the ratio 
and the portions the curve principally affected depends 
Effects 
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Finally there are variations mean sea level due periodic meteorolog- 
forces. Seasonal changes prevailing winds have measurable effect 


q 


Tidal Basin Tidal Characteristics 


December, 1959 


monthly mean sea level. and Ssa are two components designed have 
one cycle and two cycles respectively each year. The annual sea level curve 


can approximated combination these two components which have 
effect the times tides because they change slowly. Because 
slow rate angular change for these components, approximate results 
obtainable disregarding whether the epochs are Kappa, Gor The epoch Coa 

divided denotes the maximum point the annual variation sea 


level months after the vernal equinox (about March 21). Thus epoch 

represents maximum about May The epoch Ssa divided 

obtain semi-annual peak months after the vernal equinox. 

peak follows six months later. Thus Ssa epoch 140° represents maxi- 

mum sea level peaks about June and December and Ssa are 

determined from number years observation. The published values 
some places have been based one year’s observations; when this the 
case, they may unreliable. Bos 
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HARMONIC CONSTANTS 
SUMMARY 


Some the more important techniques can demonstrated using the har- 
monic constants listed below and comparing the relationships with curves 
Figs. and taken from the Coast and Geodetic Survey East Coast and West 
Coast Tide Tables. The symbols the top these figures represent from 
left right, apogee, last quarter, equator, and new moon. 


Principal Diurnal and Semidiurnal 
Harmonic Constants for Various Places 


feet deg. feet deg. feet deg. feet deg. feet deg. 


Boston -46 206 187 4.48 111 -74 .96 


Seattle 2.73 278 1.50 256 3.53 .84 341 


Anchorage 2.22 340 1.17 324 11.56 106 3.16 149 1.97 
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360° have been added obtain positive value). Going back Fig. 


have value between180° and 270°, but considerably closer the lower 
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AMPTON ROAD 
FIGURE 


Pensacola have completely different tidal regime. Although and 


are small, the semidaily components are relatively negligible. This re- 


Seattle both daily and semidaily components are large and have 
low water inequality and sequence lower low higher high higher low 


sults small diurnal tide that becomes negligible when the moon the 


equator because and oppose each other that time. 
value. Thus expect high water inequality substantially smaller than the 


lower high. Looking the curve Fig. find the above except durin 


substantial diurnal inequality. Combining epochs, Mg° 
transitional period when the moon near the equator. 


gean range than the mean spring range. 


spri 
Ft. 
| 


HARMONIC CONSTANTS 


Finally Anchorage the outstanding characteristic the large amplitudes 
the semidaily components. follows that there large semidaily range 
with relatively small inequality and, more typically than Boston, larger 
spring ranges than perigean ranges. Since 162°, the in- 
equality, although small, larger the low waters than the high waters. 

Fig. bears this out. 


The above are illustrated examples how data typical shapes the 


tide curves can readily estimated from the harmonic constants listed. 
Other deductions shape curve can made from the harmonic constants 


and Mg, and seasonal variation sea level from the components 
and Ssa. 
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INTERIM CONSIDERATION THE COLUMBIA RIVER 
Closure John Lockett 


JOHN LOCKETT,! author wishes thank Mr. 
Hickson for his discussion. During preparation the paper, Mr. Hickson was 
furnished copy the final draft for his information and took that 
tunity furnish his comments thereon. His published comments, although 
they are substantial disagreement with the views expressed the author, 
are welcomed since honest, constructive criticism believed one the 
keystones human progress. 

The paper points out that conditions the Columbia River entrance 
been undergoing and are still continuing undergo process change under 
the influence completed corrective and training works. Completion the 
original South Jetty 1895 was followed extension that structure 
1913, construction the North Jetty 1917, and construction Jetty “A” 
1939, all which were provided under the and 40-foot entrance projects. 
This, followed extensive dredging under the 48-foot project adopted 1954, 
has subjected the entrance area various degrees control, resulting 
continuously changing conditions within the entrance river and ocean forces 
attempt reach state equilibrium. Hence, within the time permanent 
improvement, not believed logical select any one short period its 
history representing optimum stage development. is, therefore, 
believed that may never possible restore entrance conditions those 
prevailing during the period 1916-32, mentioned goal achieved. In- 
stead, believed that must work with conditions they exist today and 
then determine how these conditions may further improved for tomorrow. 

Without casting discredit the engineers responsible for the design and 
construction existing Columbia River entrance structures, who certainly 
performed miracles accomplishment considering the prevailing conditions 
and the limited resources available, must admitted that practically all 
this work was, necessity, planned, designed, and constructed largely the 
basis judgment considerations. Some this work has been extremely suc- 
cessful while the true effectiveness other work still doubt. According- 
ly, was considered advisable connection with the adoption the 48-foot 
entrance project, which includes provisions for construction spur Jetty 
that further work, including dredging, based known engineering fact 
rather than judgment considerations governed the past. view the 
rapid advances made the field tidal hydraulics during the past years 
elsewhere the nation, followed that new approach the problem the 


Proc. Paper 1902, January, 1959, John Lockett. 


Head, Special Projects Investigation Section, Planning and Reports Branch, 
Eng. Div., Engr. Div., North Pacific, Portland, Ore. 


Page 
at 
S 
: 
7 


Columbia River entrance was indicated and that this new approach should in- 
clude searching engineering study all possible phenomena which might 


have important bearing the problem. Mr. Hickson’s position this re- cur 
gard, however, may summarized the following extracts from his initial itr 
comments the final draft the paper: ary 
appears are making very complicated problem out one 
which relatively simple. simply problem river hydraulics 

which the principles hydraulics apply. differs from inland river 
control only the extent that winds, waves, tidal currents, density 
currents, and littoral currents distract our attention from the basic 
consideration, which that the river outflow will produce and maintain 

satisfactory channel properly directed alignment and width 
between controlling bank constructed works. This has very success- ties 
fully been done the upriver ship channel the Columbia and other 
solution made more difficult attempts make detailed analysis the 
all the forces which are work and what part each may have pro- 
This concept raises question how may properly direct the river gre 
flow produce and maintain satisfactory entrance channel when know cur 
neither the full nature nor magnitude the forces controlling the regimen 
the estuary, let alone the effects these forces any further regulatory vati 
structure which may devised. 
The author, general, does not believe that any strong exception could diti 
taken the discussion provided was based conditions prevailing prior 
improvement the Columbia River entrance when the forces nature were and 
such prevent saline intrusion into the estuary. The subsequent con- 
struction corrective works and dredging, has, however, upset the balance 


point where they are believed represent the major dominant force control- 
ling the regimen the estuary today. This logical sequence events 
which has been confirmed histories estuary improvement elsewhere. 
There are, however, several points brought out the discusser which re- 
late quite fundamentally the principles tidal hydraulics and should, 
fore, considered further detail. First all, not believed that 
design any estuary structure can ignore saline influences represented 
density currents. These currents are created differences specific 
gravity the fresh river water and the salt water the sea. The magnitude 
these currents difficult fully describe but may visualized some 
extent consider the hydrostatic conditions created the meeting the 
fresh river water weighing about 62.4 pounds per cubic foot with the salt water 
the sea weighing about 64.0 pounds per cubic foot. The difference spe- 
cific gravity, reflected the difference weight the sea and fresh water, 


} 
natural forces and has allowed the penetration saline influences 


1.6 pounds per foot depth, creates hydrostatic head upstream 
rection with force about pounds per square foot the 50-foot depth 
level. This equivalent total resultant force over 5,000,000 pounds Fig. 
when applied through the authorized one-half mile channel width the 
Columbia River entrance. The true magnitude this force becomes Gua 


what greater than 5,000,000 pounds when consider that the full width be- 
tween the Columbia jetties least two miles. This force, evidenced ther 
density currents, acts throughout the region saline penetration within the 


3 
% 
4-2 
rs 
‘ 


DISCUSSION 


estuary and kept balance counteracting forces friction along the 
estuary bottom and inertial energy flowing fresh water. Until this density 


in- 


re- current force finally dissipated the upstream limit saline penetration, 
remains factor deserving detailed consideration the planning estu- 
ary structures. Therefore, the author cannot agree with the concept that 
should not allow saline influences distract our attention the consideration 
Considerable emphasis has been placed the discussion the interpreted 
results the 1932 current survey which was undertaken along single range 
across the estuary about five miles above the outer ends the Columbia 
jetties. Equipment employed measure current velocities different depths 
each the five measurement stations along this range consisted the 
ordinary Price current meter suspended from the metering vessels. Facili- 
ties determine measure the direction subsurface currents were, 
course, not available that time and, consequence, the actual direction 
subsurface currents was largely matter personal judgment. Further, 
the survey did not include concurrent measurements salinity and tempera- 
ture along this range and, hence, was not possible confirm the judgment 
decisions current direction with measured indications salinity. re- 
grettable that substantial part the detailed field records made during this 
current survey has been lost over the years and, accordingly, its value has 
been reduced circumstances beyond control. Even though all these obser- 
vations were available today, not believed, due changes occurring with- 
the interim period, that they could accepted representative con- 
ditions presently prevailing the estuary. 
this connection, the discusser states that, general, the same littoral 
were and density currents prevailed the entrance during the period 1916 1932 


prevail today. Although the magnitude the littoral effect still un- 


measured and, therefore, the accuracy this part the statement can neither 
confirmed nor denied, there logic the belief that the density current 
pattern substantially stronger today than was during 1916-1932. Exami- 
nation condition surveys undertaken during that period reveal that control- 
re. ling depths the entrance, except between 1925 and 1927, were less than 


re- feet and 1916 were only feet. 1925 and 1927, controlling depths were 
and feet, respectively. recognized principle tidal hydraulics 
that when depths estuarine channels are increased for navigation other 
nted Purposes, the effective density head the entrance increased and density 
action within the estuary becomes more pronounced. Since 1932, the 
mitude controlling depth between the jetties has been increased least feet and 
some is, therefore, logical expect greater density current action now than dur- 
the ing the referenced period. Even though the channel section were held constant 
the construction dikes contractive works, such the proposed Jetty 
deepening the entrance channel must increase density current effects 
water, the saline region the estuary. 
di- The discusser also states that salinity measurements through tidal cycle 
show that most the salt water washed out each ebb tide. refers 
his discussion showing graph form measurements specific 

taken opposite Fort Stevens, Oregon, near the Point Adams Coast 
ome- Guard Station, river mile 7.5, 1936 support for this conclusion. Although 
known that actual density water varies somewhat with its temperature, 
inthe time. Therefore, the absence temperature data, doubted that 


: 

§ 


these observations can considered conclusive. any event, difficult 
reconcile these observations with those made practically the same lo- 
cation 1955 disinterested agency reported Table the paper 
under discussion. noted that, forwarding his report the 1936 obser- 
vations the Committee Tidal Hydraulics, Corps Engineers, Mr. Hick- 
son, his letter memorandum April 10, 1953, stated that re- 
alized that the basic data secured were not comprehensive enough admit 
drawing firm conclusions”. Some six years later, however, these same obser-| 
vations are presented firm data warranting the conclusion that ebb forces 
remove most the saline influences from the entire estuary. work has 
been undertaken since 1953 which might lessen the magnitude saline influ- 
ences the Columbia River estuary and there appears basis for re-| 
appraisal the value the 1936 observations. Even though these particular 
observations were such nature permit the formation firm con- 
clusion, such conclusion could only apply conditions prevailing the vi- 
cinity Fort Stevens, river mile 7.5, but could not indicative conditions 
occurring elsewhere the estuary. 

The discusser also stresses the need for early construction spur Jetty 
“B” regulatory structure designed reduce dredging and provide stable 
channel through the entrance. presents cross-sections the estuary 
the general vicinity this proposed structure and indicates thereon the ap- 
proximate flow area that would cut off the structure. the opinion 
that Jetty “B”, denying the river use the area thus cut off, would con- 
tract the available cross-sectional area the point where increased river 
velocities would solve the shoaling problem erosion Clatsop Spit shoals 
along the opposite shore. 

Examination the accompanying Fig. shows that proposed Jetty “B” 
would lie along the north shore the estuary, almost parallel the North 
Jetty, with its extreme end approximately bisecting the distance between the 
ends Jetty “A” and the North Jetty. Fig. indicates that Jetty “B” would 
constructed depths averaging less than feet below MLLW and maxi- 
mum depth feet its southernmost extremity. total cross-sectional 
area about 521,000 square feet, Jetty “B” would occupy about 160,000 square 
feet, roughly per cent the total section. Under conditions normal 
river flow, where there are saline influences reversals currents 
distract our attention, such reduction section could expected produce 
some effect the direction river flow increasing current velocities 
the point where they might erode less firm areas along the opposite shore. 

Considering, however, flow conditions they are believed exist this 
part the Columbia River estuary, the prospect solving the problem 
this relatively simple expedient becomes somewhat less than promising. 
due the fact that reversals flow, created tidal action, introduce 
factors that complicate nullify the application the principles pure 
hydraulics. Under such conditions, may seen that the confining effect 
Jetty “B” would limited that portion the time which the entire flow 
throughout all depths this restricted cul-de-sac area, the 
Columbia channel the south, would ebb direction. believed 
that this condition, however, could expected prevail for only relatively 
short periods time due bottom flood flows generated density current 
action, even during ebb tide. During that portion the time which the 
flow the cul-de-sac would flood direction, appears that Jetty 
would unable affect the critical portion Clatsop Spit lying the west. 
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The structure might, instead, create eddy conditions which could conceivably 
threaten the stability the North Jetty. Further, the main pattern flow 
this portion the estuary would continue unaltered the area immediately 

south the structure where depths extend more feet. believed 
that any density current action occurring that principal flow area would, 
therefore, unaffected Jetty “B”. Hence, there appears exist reason- 
able basis for doubt regarding the true value that structure. 


The discusser claims that construction Jetty “A” 1939 stabilized the 
| 


channel that point, caused widening the channel for 800 feet the 
south, and presently maintains depths excess feet for distance 
over mile downstream. Although the author does not disagree with evidence 
showing that the 40-foot contour opposite Jetty “A” moved southward following 
construction that structure, must pointed out that this local improve- 
ment was accompanied deteriorating conditions the northwestern portion 
Clatsop Spit, where the 40-foot contour that area advanced markedly 
toward the north and west. Also, examine condition surveys prior 
1939, will noted that depths considerably excess feet have pre- 
vailed for over mile downstream from this point ever since improvement 
the entrance was initiated. Accordingly, not clear how Jetty “A” may 
properly credited with maintaining depths which have always existed that 
area. 
regard the effectiveness Jetty “A” and the dikes along Sand Island, 
and the need for comprehensive study the problem, the Committee Tidal 
Hydraulics, its reply June 12, 1958 Mr. Hickson’s comments Com- 


mittee Technical Bulletin No. which are both listed unpublished refer- 
ences the discussion, stated: 


The second exception concerns your statement that the spur dikes 


since their completion about 1938; the Committee does not share your 
view that these structures have been effective claimed. The 
Committee feels very strongly that the recommended program proto- 
type measurements and the model study the Columbia River entrance 
should carried out the earliest practicable date. Until these 
missions have been accomplished, and the results obtained have been 
fully analyzed and evaluated, the Committee believes that the con- 
struction any so-called improvement works the problem area 
would grave mistake. The existence such wide differences 
opinion the basic reasons for this shoaling problem, among engi- 


along Sand Island and Jetty “A” have maintained excellent channel 


neers who have devoted most not all their careers the study 
tidal hydraulics and related matters, ample proof the Committee 
that this problem worthy careful and complete study and evaluation 
prior the expenditure large sums for improvements”. 


Fortunately, the plan action recommended the Committee was initiated 


Portland District with the completion 1959 prototype measurement 


program which observed current velocities and directions well salinities 


and temperatures some stations along six ranges across the estuary. 
Observations these phenomena were made five depth levels minute 
intervals each station throughout complete tidal cycles hours each 
under conditions low, normal, and high river discharge. The observed data 
are currently being evaluated and the author hopes present some the 
the program forthcoming paper. 


» 


af 
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Preliminary the conduct this measurement program, Portland 
made several tests current velocity and direction equipment offshore the 
Columbia River entrance. had been suspected that dredge spoil 
Areas and which had been used for this purpose for many years, 
its way back into the entrance channel, tests were made the vicinity 
particular spoil disposal areas. These tests revealed that such was quite 
likely the case and, beginning 1958, Areas and were abandoned and 


practically all spoil was disposed Area depths 120 feet more. 
This change disposal practices considered the author responsi- 
ble for the securing 1958 and 1959 improved entrance depths the end 
each dredging season substantially reduced dredging costs. 
that through this single expedient, the savings dredging costs thus obtained 

will repay few years, the entire cost the Columbia River tidal hydraulic 
investigation. 


The author recognizes that quite probable that final evaluation the 
prototype measurement program may not completely confirm all his views 
garding the behavior currents and other phenomena the Columbia River 
entrance. That little importance. The important thing is, course, the 
fact that positive action now being taken learn specific details regarding 
the forces controlling the regimen the Columbia River estuary that logi- 
cal and factual planning can precede actual design and construction any 
further regulating structure structures that may found necessary. 

Adoption modified spoil disposal practices the Columbia River en- 
trance, based modern technological knowledge the field tidal hy- 
draulics, has already been rewarding and the continued application this 
knowledge the over-all problem holds the promise still greater rewards 
the future. 
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HYDRAULIC ANALYSIS SURGE TANKS DIGITAL 


Discussion John Herak and Salvador Rodriguez 


JOHN HERAK! and SALVADOR ASCE.—The author 
has made important and significant contribution the literature hy- 
draulic transients. The writers know other contribution opening such 
large avenues refinements the treatment transients and eliminating 
the usually tedious labor making arithmetic integrations. Two major points 
impress the writers. The first the close check given the manual compu- 
tations for all large maneuvers. This also emphasizes again the inherent 
safety resulting from the use orifices equivalent restrictions. The 
second the opportunity for making refinements the solutions, thus pro- 
ducing mathematical results for many minor factors that heretofore were 
necessarily intelligent guesses, unless undue amount labor were 
expended. the other hand this requires closer look the assumptions 
and coefficients used for these minor effects, order for them not mis- 
leading. hoped that Yu-Tek Li’s researches, referred later, can 
published soon, since they indicate considerable variation from some the 
results, such those Vogel and McNown, previously available. 

The following evaluation the differences between manual and ma- 
chine results the surge tank stability study. 

The evaluation tee losses, the arithmetic integration study, were 
based upon values obtained using Vogel’s data for shape (rounded edges) 


tee connections, and branch main diameter ratio 0.58. Tee loss coef- 


ficients prepared the Omaha District, C.E., (Table from Vogel’s data 
were based branch main diameter ratio 0.667 and zero discharge 
ratio have been set equal zero, thereby modifying Vogel’s results the 
range zero 0.1 discharge ratios. The effects these differences be- 
tween manual and machine coefficients are noted below. 

The arithmetic integration method investigating the stability the Oahe 
System was subject the unit relatively small instantaneous load change 
and determine whether the amplitudes the tank levels with respect the 
steady level were damped magnified with increasing time. stated 
the author, manual stability computations proceeded directly from point 
point (Fig. 11) and once entry onto the curve was made, constant 
power was assumed delivered the turbine full gate position. 

This transient produced oscillations and along the constant curve 
78000 without any time reaching full gate. The arithmetic integration 
neglected the effects governor dead time, friction losses, inertia, and 


Proc. Paper 1996, April, 1959, Nicholas Barbarossa. 
Princ. Engr., Sverdrup Parcel Eng. Co., San Francisco, Calif. 
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elastic energy the hydraulic system downstream the surge tanks; and the 
energy storage and inertia the rotating mass and electrical system. The 
computer studies included only partial corrections for some these factors 
and did not consider the others. However, step the right direction was 
made including the partial corrections. 

partial arithmetic integration study was made Test No. Fig. 11, 
ing time interval 0.1 second and the results are shown Fig. 
The transient path full gate position similar that obtained the 
computer, the variation apparently being due differences values 
loss coefficients. both cases, the transient attains full gate position before 
reaching the curve. believe the computer should have been pro- 


grammed permit the transient follow the full gate curve, indicated 
Manual study intersection with with some overshoot due gover- 
nor dead time and system inertia. (Dead time considered equal 0.05 second 
for Manual Study this point the governor takes over and begins 
the turbine gates since the turbine output exceeds the required horsepower. 
The path the transient will approximate spiral and will converge the 
new steady state point, the unit operation stable. suggested that the 
computer programmed outlined above and allowed traverse through 
several cycles obtain better measure both surge tank and unit operating 
stability. would interest compare the decrement (Jaeger), 
measure the damping the hydraulic system, obtained machine 
with that obtained the manual method. 

The transient curves Fig. all indicate that the gates begin close 
before the new steady state power attained. This not compatible 
with governor action. 

previously noted, differences tee loss coefficients existed between 
manual and machine computations. recent unpublished paper entitled 
Loss Tee Section and Manifold”, Yu-Tek Li, February 1959, former} 
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the member our hydraulic section connected with the Oahe Reservoir Project, 


presents values head losses tee sections, head losses through orifices 
the case where the plate orifice installed the lateral pipe, including 


values for different shapes orifice plates, and for different locations the 
orifice plate the lateral pipe; together with determination head losses 
manifold system. Attention being called the paper, not with the in- 


tent defending either set tee loss coefficients used these studies, but 

With the thought that while the writers agree with the author’s ideas the 

solution transient problems computer, additional basic work still re- 
evaluating items such head losses tees before “exact solutions” 


can made with computers. 
conclusion, the writers wish commend the author his excellent 
paper and look forward the publication his paper “Proposal for Compre- 
hensive Study Digital Computer Hydraulic and Governing Transients 
Affecting Design and Operation Power Plant System”. 
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RESISTANCE EXPERIMENTS TRIANGULAR 
Discussion van Malde 


van paper deserves study the experiments described 

are unique regarding the sizes the flume used, whilst the accuracy with 
which they were carried out make the results most valuable. The paper af- 
firms that open channel cannot calculated with the well-known formulas 
Prandtl (which were found from experiments with circular pipes) without 
making error. According their interpretation the data the authors 
state that the average error the mean velocity when using these formulas 
would amount 3,5% for the observed triangular channel being smooth 
flume and 2,2% the case the channel being rough flume; these numbers 
refer tranquil, non-critical flow. The authors conclude that both cases 
the logarithmic formulas Prandtl have advantage over the Manning for- 
mula and that the second case the formula Morris for isolated-roughness 
flow has advantage over Manning’s Prandtl’s formulas either. 

studying the data for the hydraulic rough flume (runs 49-101) appeared 
that the authors put together all values (the relation between observed 
and for smooth pipe) for certain roughness, irrespective slope (S) and 
depth Kirschmer, who far the writer knows first pointed the 
resemblance resistance hydraulic rough open channel hydraulic 
smooth circular pipe mentions one his papers(2 the possibility that 
influenced This concludes from some the experiments 
Bazin, taken rectangular flume with battens fastened the walls the 
direction the flow; would increase with increasing Theoretically there 
explanation for this influence, but the theory for the uniform steady 
flow open channels still doesn’t take into consideration all the phenomena 
(such secondary currents) there reason exclude the possibility 
this influence. For the same reason may possible that the depth flow, 
Yn, influences (though this seems unlikely directly related the 
Reynolds number, which used compute 

Fig. shows (calculated with slide-rule accuracy) function for 
different values From the points for runs 84-92 and 93-101 may 
concluded that for these cases independent Yn; the other series 
points show trend for decrease with increasing depth, but the points 
generally scatter too much take this for granted. the writer’s opinion 
the points definitely show that increases with increasing slope; this the 
Same result Kirschmer did find. 

smooth tranquil flow the authors found the indeed astonishing 
result that 0.93. For the runs concerned (1-16; 102-108) was again 


Proc. Paper 2018, May, 1959, Ralph Powell and Chesley Posey. 
Civ. Engr., River Research Section, Directorate Upper Rivers, 
Rijkswaterstaat, The Netherlands. 
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plotted against for different (Fig. 2). The figure reveals that the low 
average mainly determined low values whilst there seems 
tendency for increase with increasing slope. doesn’t seem justified 
state anything more definite, but the trend the same the one for the hy- 
draulically rough flume. 

The values found here make possible find two equations 
application Morris’ theory(3) considered suitable: one for 8.55 
and one for 19.3 The relevancy Morris’ formu- 
las for practical computations doubted the writer they seem too com- 
plicated; empirical relation will suffice for “isolated-roughness flow” 
and “skimming flow”. 

The results this analysis indicate that the conclusion the authors re- 
garding the equal relevancy the formulas Prandtl, Manning and Morris 
for rough triangular channel may premature. The writer had oppor- 
tunity deal with this question detail; maybe definite answer can only 
given after new series tests will have been carried out even the ex- 
cellent data given the authors show for some runs (57-64) considerable 
scattering points. This question equal relevancy formulas may seem 
matter pure theoretics without much practical meaning for the flume 
studied, the mean discrepancies the mean velocity all three formulas 
only amounts 2,3%, whilst, according the authors, the probable per- 
centage error determining the discharge for some the runs may amount 
the same percentage. When studying river problems however the influence 
may importance and for this reason attention drawn some 
conclusions drawn from further analysis the data given the authors. 


REFERENCES 
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Morris: New Concept Flow Rough Conduits. Proc. ASCE, 
Vol. 80, nr. 390, January 1954. 
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RESISTANCE PROPERTIES SEDIMENT-LADEN 
Discussions Tsung-Lien Chou and Simons and Richardson 


TSUNG- LIEN ASCE.—The authors are highly compliment- 
for taking important but long neglected problem into the laboratory and 
the library with special technical dexterity and mathematical insight elimi- 
nating the interwoven elements and bringing the required factors focus. 
The results are neatly tabulated and the conclusions drawn are concise and 
precise. Since the problem not only important for practical hydraulic engi- 
neers dealing with river and canal works, but also closely related with 
the more theoretical and broader side boundary layer and turbulence 
general, the writer would like bring out some points which may help clarify 
the issue. 

Table shown that each set observation, all important hy- 
draulic characteristics are identical except the slope and the sediment dis- 
charge concentration With these two, the authors attribute direct cause 
effect correlation friction factor repeated Table the 
figures are consistent and impressive, though not quite proportional. What 
this amounts that the structure this part the problem depends the 
accuracy and reliability the measurement which always delicate and 
complicated. this junction, may mentioned that empirical formulas 
used fit the observed values and they are supposed take care all 
losses which will, modern hydrodynamical terminology, include losses clue 
skin friction, form resistance and turbulence the pure fluid side and loss- 
knocking out solid particles from embedment, bringing them into sus- 
pension, and accelerating them the same velocity the dispersion medium, 
and host other losses which are yet unknown. widely observed, the 
velocity varies with space and time, even uniform conduit under steady 
flow. With the empirical relation and indicated above, synchranically 
together with varies the slope which may include many ups and downs, 
even negative values backflow pockets. natural streams, only mean 
values are taken for both. Therefore, real sense, both and are statisti- 
cal nature and nothing more than mathematical fictions. Similar are the 
friction factors, Darcy-Weisbach for Manning Furthermore, these factors 
are not universal constants shown the following tabulation. Thus, each 
friction factor has its own structure and the relations among them are not 
homogeneous. The fact that they can used calculate flow their special 
forms does not necessarily mean they can represent certain factors 


Proc. Paper 2020, May, 1959, Vito Vanoni and George Nomicos. 
Hydr. Engr., Clinton Bogert Engineers, New York, 


to 

> 

ume 

r- 

nt 

ice 
me 

iN 


108 


Manning 

that 

Original Forms vezcyYrs It s 

1.36 whethe 

Energy Slope s 1,667 ¢ 71.167 2.203 373 


have 
for 


For equal s, 


Ratios of Factors 


0.091 0.0381 
1.372r v 


suctio 
modern hydrodynamics. Now Table held count for the variation ™°re 
minute vortices and others. This apparently outside the function assigned The 


applying the empirical formula flume like this, should cover only 

that portion with established normal flow, and with the entrance and tailwater 

regions deducted and the cross-section observation would set the direct 


bed the differences water levels the ends the flume may bring 
drop-down back water effect, which would not represent the true value 
solidifying the flume bed after running with sediment, the authors tried factor 
eliminate all factors due bed configuration order emphasize the 
influence sediment hydraulic resistance. may worthwhile recall Fir st, 
that the turbulence produced movable bed not the same 


geometrical center the portion. Any average value from tilting flume fact 


bed which fixed, even though the configurations are the same. 
bed, turbulence evened out picking solid particles points strong 
intensity and dumping some its load points weak intensity, the stage. 
distribution turbulence contours more less conformal the bed to- 
pography. Contrarily, fixed bed, irregular configuration tends amplify This 


the turbulence points strong intensity and therefore greater turbulence is} the 
expected fixed bed. Thus some part the calculated value might offset this 
this cause. the 

The velocity profile shown Fig. are apparently not taken from the ex- with 
periments this flume (10-1/2 inch wide). the depth, the slope Provi 
width are identical for the two profiles, the one with higher velocities will loess, 
higher mean velocity and thereby bigger discharge. that case, more power 
required the and there hardly any strick comparison the 
two. 

The authors maintain that the suspended sediment reduces hydraulic basin. 
ness its damping effect turbulence the following reasoning: “To 
keep sediment suspension, i.e. prevent from settling due gravi- 
tational force, work must come from the vertical component turbulence 
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fluctuations and must result damping turbulent motion” and “this means 
that the momentum transfer coefficient also decreased thus allowing the 
velocity and velocity gradient increase.” 

should remembered that turbulence created the expense 
energy stock the flow. That say, soon the turbulence gener- 
ated, certain amount energy lost from the flowing stream, matter 
whether energy used keeping sediment suspension supplying 
Vortex motion minute masses and finally dissipated into heat. keep 
sediment continuously suspension, uninterrupted supply the necessary 
energy through turbulence required. However this damping action may 
have the effect preventing the eddies spread out into the entire stream 
localizing them the lower strata. Finally stream where the saturated 
suspension tends dump its load, there chance releasing some energy 
for accelerating the flow. These conditions are limited certain part 
certain period flood where suspended load dumped. 

Furthermore, the transportation processes solid particles are different 
for different effective sizes and distributions. Only fine particles less than 
0.2 can thrown into suspension and move with turbulence gener- 
ated purely hydraulic resistance. the other hand, bigger particles cre- 
ate extra turbulence their projections into flowing stream and produce 
suction effect helping them into sliding, rolling and saltation. This demands 
more energy consumption and thereby increases the hydraulic resistance. 

The above conditions are well observed pumping sands suspension and 
suspension tests. 

addition the direct experimentation laboratory, the authors collect 
some interesting information field data from national streams. Here the 
direct cause-effect correlation and further overshadowed addition- 
factors (a) source solid supply, varying with the topographical and 
geological conditions the river valley, (b) hydraulic properties the 
stream, aggrading degrading. The mutual action and reaction these 
factors can easily produce contradictory effects. 

general, rising stage has higher average velocities for two reasons. 
First, rising stage, the energy slope drop-down curve with increasing 
velocities towards each downstream section. Thus the velocity accelerated. 
Secondly, the potential head converted into kinetic head with smaller loss. 

most cases, the amount and the intensity sediment load vary with the 
stage. freshet originated from upstream watershed and gains 
strength flowing downstream, will pick silt its optimum capacity. 
This happens very often the rising stage. After passing over the peak flow, 
the river tends dump its suspended load. However, there may exceptions 
this rule dealing with difference sources solid supply. For instance, 
the Yellow River Tung-Kuang, the intensity suspended sediment varies 


with the source flood. comes from the north bank tributaries (Shansi 


Province) where steeper grades river beds with thicker eolian deposits 
loess, the River has tendency dump its load even rising period. The 


true, the flood started from the south side. Similar situation 
observed Chingling Chi the Yangtze River where the discharge from 


the south clarified passing the Tung-Ting Lake which acts settling 


basin. 


The above few examples just illustrate the complicated relations various 
factors natural streams and the cause-effect correlations can easily 
reversed, unless comprehensive assessments are made for some period. 
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the theoretical side, the sediment can reduce hydraulic resistance rela 
the following ways: 
abov 

(1) Smoothing the surface roughness scouring off the mound and pres 
the depressing the 

(2) Evening off the aquatic growths flum 


(3) Exposing fresh surfaces which are more pliable 
(4) Compacting particles selective deposition and distribution ness 
(5) Keeping down the vertical component turbulence (as the authors’ 
q 


limi 

main item). Again, every factor must properly appraised. size 

The problem suspended load hydraulic resistance complicated but the 

very interesting. The authors have opened new vista for both practical 


and theoretical investigation. Perhaps, another way attack may the di- 
rect measurements the turbulence itself for the similar condition turbid 
and clear flows. 


authors are commended for their excellent explanation many the 
important aspects the mechanics flow alluvial channels. Based upon 
similar studies conducted both large and small flumes the writers, the 
equipment used and few the concepts presented the authors are elabo- 
rated on. 

One the problems associated with study this type the difficulty 
evaluating the effect flume size data and observed flow phenomenon. 
The writers found that the results obtained with given bed material using 
large recirculating flume 150 feet long, feet wide, and feet deep, could not 
reproduced smaller recirculating flume feet long, feet wide and 
feet deep. For example, the forms bed roughness observed the order 
increasing shear the small flume were: 


Plane bed before beginning motion 

Plane bed after beginning motion 

Dunes 

Transition from dunes plane bed and standing waves 

Standing sinusoidal sand and water waves which were phase 
Antidunes 


Using the same bed material the large flume the forms bed roughness 
order increasing shear were: 


Plane bed before beginning motion 
Ripples 

Ripples superposed dunes 

Dune 

Transition from dunes plane bed and standing waves 
Standing sinusoidal sand and water waves which were phase 
Antidunes. 


The large flume produced the forms bed roughness observed the 
whereas, the small flume did not. was also noted that the height dunes 


Assistant Prof., Civ. Eng. Dept., Colorado State Univ., Fort Collins, 
Hydr. Engr., Civ. Eng. Dept., U.S.G.S., Fort Collins, Colo. 
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related the depth sand bed. depth bed material decreased the 
amplitude the dunes decreased. These differences results, over and 


above those which can compensated for correcting for side wall effect, 


present serious problem those who are trying utilize data collected 
the various investigators, particularly that data gathered very small 
flumes. 

mention made the authors the fact that the form bed rough- 

ness related size and gradation bed material. This may due the 
limited range conditions which they investigated. illustrate using sand 

sizes ranging from 0.2 0.5 the writers have observed that 

the amplitude the dunes independent size bed material but the spac- 


cal ing and shape the dunes are related size bed material. The spacing 
the dunes increases with decreasing size and the angle the fore plane the 
irbid dune makes with the horizontal decreases with decreasing size bed materi- 


al. The net result that resistance flow smaller with the fine sand than 
with the coarser sand much per cent. This discussion perhaps 
raises the question bed form terminology. avoid confusion, the major 
forms bed roughness which have been observed alluvial channels the 
the writers are illustrated Fig. 1.* From this figure, apparent that the 
upon form bed roughness described the authors, according this classifi- 
the cation, ripples. 
Referring the difference the bed friction factor for sediment-laden 
flow over loose sand bed and clear water flow over fixed sand bed, the 


same configuration, doubtful that the total difference can explained 

order 
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the presence the suspended sediment. The writers have observed that the 
scale and intensity turbulence are entirely different rigid bed with 
bed material transport than loose sand bed having the order ppm 
bed material transport, none suspension, and other flow factors the 
same. 

The effect sediment load resistance flow significant illustrat- 
ed. However, the authors conclusion that all the computed decrease re- 
sistance flow can attributed the damping effect the suspended sedi- 
ment the turbulence not concurred in. The presence fine sediment 
causes effects fluid properties which change the entire fluid bed material 
interrelationship. These factors appear more importance than the 
dampening effect the suspended sediment. 

The presence very fine sediment (clays) the water also influences re- 
sistance flow alluvial channels. Based upon recent studies the 
concentrations fine sediment the order 40,000 ppm can reduce re- 


sistance flow much per cent the tranquil flow regime. the 
under rapid flow conditions, the presence fine sediment can increase re- 
sistance flow. the 
The validity the page 104 which explains the existence 
smaller friction factor for the same discharge rising stage than the 


falling stage, where the decrease friction factor presumed result 
readjustment the bed due the larger sediment load questioned. 
not more probable that the increase load and the reduction friction factor 


are both related the fact that with the increase shear with stage the bed 
form changes and that the change bed form will lag the change stage? 

better cope with the problems associated with fully developing and the 
izing resources, additional work this caliber, preferably done 
larger flumes, should initiated. Only making studies this type can achie 
hope explain the quantitative effect size bed material, gradation 
material, temperature, geometry channel, depth sand bed, fine 
load, and rate change discharge with time total bed material transport, 


the forms bed roughness which develop, and the related resistance flow. 
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HYDRAULIC MODELS THE ST. LAWRENCE POWER 


GEORGE BERRY,! ASCE.—The St. Lawrence Power Project was 


unusual and challenging many respects. Although the St. Lawrence River 
normally reaches its minimum annual discharge during January February, 
the average January discharge for years record 221,000 cfs and the 
average February discharge 218,000 cfs. Comparing these averages with 
the average annual discharge 246,000 cfs apparent that low flow 
season for diversion closure does not really exist. divert such large 
flows required that the construction performed several stages. (1) 
each construction stage the depth water the navigation channels and 


canals the project area had maintained. The water surface elevation 
Lake Ontario could not altered from its natural condition and the dis- 

charge downstream power plants could not impaired. 

The hydraulic model studies described Mr. Bryce focused attention 

the salient features possible alternate solutions and permitted the hydraulic 


engineers experiment with these variations until the optimum design was 
achieved. 
This series model studies outstanding example the technique 
complementing the general results small scale model investigation 
rt, the critical details larger scale models. The best location for the 
structures and the location, size and shape the channels required navi- 


These were the scale 1:500 horizontally and 1:100 vertically which 
quite large for models such large river. Details the stilling basins and 
gates were then studied 1:50 scale partial models. Finally comprehensive 
models Long Sault and Iroquois with enough upstream and downstream river 
channel simulate prototype approach and discharge conditions were built 
and tested determine the flow pattern and hydraulic performance for each 
construction stage. 
the energy dissipators Long Sault and Iroquois were developed 
from long series tests for each. Several forms baffles and sills were 
tried. The final form which had baffles and solid end sill was the results 
careful experimentation obtain the lowest possible velocity the 
channels below the dams. The ineffectiveness baffles these spillways 
Which have high entrance velocities agrees with experiments conducted few 
years ago M.LT. the effects Baffle Piers The ex- 
cellent performance the Long Sault stilling basin can judged from the 
fact that, although the natural tailwater depth was insufficient form 


gation and ice cover considerations were determined the river models. 


é 
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hydraulic jump, jump was forced form within the basin for the three 
different diversion conditions (a) diversion through low sluices, (b) di- pape 


versions through ports left the dam and (c) flow over the 


The channel downstream from the dam was rock which the natural pre- invol 
project velocity was the same range those achieved with the recommend-§ 
spillway design. Similarly satisfactory energy dissipation was achieved 
Iroquois Dam. 


The remarkable reduction the hydraulic downpull effected reshaping featu 
the bottom the Long Sault port gates illustrates how important this factor 
designing gates and their hoisting equipment. Unbalanced hydraulic ately 
gates this type depend primarily the local conversion pressure 
head velocity head the restriction beneath the gate. The shape the hydr: 
bottom the gate and the approach conditions determine the magnitude the 
pressure reduction and the area over which acts. Size gate, depth sub- 


mergence, and location skin plate and seals determine the magnitude the 
hydrostatic force which not counterbalanced vertical force the bottom 
the gate. The tests described cover variety gate conditions with up- north 
stream and downstream skin plates and with various depths submergence. the 
believed that structural and hydraulic engineers would find the data very Eiser 
valuable the author would show the physical details and dimensions the 
gates, the pertinent boundary conditions and the head the gate along with condu 
uplift and/or downpull various gate openings. Available data correlating velop 
these factors are scarce. quenc 
impossible accurately isolate the savings construction cost at- strea 
tributable directly the model studies. The substitution dry excavation 
south Gallop Island for wet excavation along the north channel resulted this 
saving about $1.00 per cubic yard almost million cubic yards exca- the 
vation. This alone repair manifold the cost the model studies. equal chann 
importance was the fact that these accurate model tests, the design was 
ties for both navigation and ice packing criteria could met with confidence dike 
that very little margin for unknown unexpected conditions was needed. 
The definite detailed construction procedure which was developed and reach 
the specifications and the model demonstrations, which were some dey 
made for the bidders show the construction sequence and hydraulic con- cept 
ditions which would encountered, undoubtedly resulted substantial 
cost. distor 
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paper portrays classic example the utilization hydraulic models 
fundamental part the design and construction major engineering project 
involving wide variety complex interrelated features. Wide experience 
has established the fact that model studies are indispensable adjunct 
analytical procedures the development optimum designs and construction 
sequences for major hydraulic channels and structures involving complex 
features and large expenditures funds, and that model studies almost in- 
variably produce major design improvements costs that are disproportion- 
ately low relation project costs. For many hydraulic engineering 
problems, including some the features the St. Lawrence project, 
hydraulic model studies constitute the only reliable approach satisfactory 
solutions. 

noted from Mr. Bryce’s paper that his lowermost river model (Model 
No. 9), utilized for studies the powerhouse tailrace, terminated just below 
the head Cornwall Island and included only the upper ends the channels 
north and south that island. mentions that portion the excavation 
the south channel, which contains the navigation channel approaching the 
Eisenhower Lock was accomplished the United States Seaway Development 
Corporation. might interest note that two river model studies were 
conducted the Army Engineer Waterways Experiment Station de- 
velop optimum excavation and spoil-placement plans and construction se- 
quences for the navigation channel the river from Barnhart Island down- 
stream Lake St. Francis. 

interesting scale-distortion problem connection with the modeling 
this reach arose from the fact that per cent the river flow turns sharply 
the south around the head Cornwall Island and enters the navigation 
channel below the lock almost normal the east-west navigation course. 
was therefore necessary, among other things, design curved deflecting 
dike realign this southward flow with the east-west navigation channel 
eliminate intolerable cross currents the latter. One model the overall 
reach, with linear scales 1:300 horizontally and 1:100 vertically, was used 
develop plans throughout the Barnhart Island Lake St. Francis reach, ex- 
cept the vicinity the above-mentioned sharp changes flow direction 
near the head Cornwall Island. Since was considered that model scale 
distortion would introduce significant errors the latter area owing the 
sharply changing flow directions, this immediate area was duplicated an- 
other model with undistorted linear scale 1:100. was expected, com- 
parisons detailed data from both models demonstrated that the distorted 
model produced erroneous distributions velocities the vicinity the 


Sharp direction changes while data other areas checked well both models. 


Accordingly, the two models were used together very effectively develop 
plan details for the overall reach. 

The problem energy dissipation downstream from the St. Lawrence 
structures unique that large discharges must passed during the con- 
Struction stage. many other rivers, seasonal fluctuations are such that 
critical stages construction can accomplished when river flows are low. 
Because insufficient tailwater was necessary force hydraulic jump 
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for diversion flows and for final operating conditions. the case the Long 
Sault Dam, this appears have been accomplished use 12.5-ft-high 
end sill. However, question arises the necessity for 40-ft-wide sill. 
lesser width sill and longer apron elevation 140 might improve 
Provision steeper slope the toe the overflow section sloping the 
diversion ports are other alternatives for improvement energy dissipation. 

The energy dissipator for the Iroquois Dam appears designed only for 
diversion flows. was stated that only ten sluices (or gate bays) were pro- 
tected, which raises the question the treatment the bed the stream 
immediately downstream from the other sluices. economies were affect- 
the omission any protection below sluices, the author com- 
mended. Many times engineers are reluctant omit any part the energy 
dissipator, even where conditions are such that can done safely. 

The effect gate operation energy dissipation mentioned only brief- 
Mr. Bryce. However, various combinations gates and openings can 
have important bearing the strength and magnitude river currents, 
and the design energy dissipator. many river dams, gate oper- 
ation secure satisfactory energy dissipation one the most important 
phases model study. 

Another important phase the St. Lawrence work reported upon Mr. 
Bryce concerns the hydraulic loading gates. This topic could the subject 
entire paper. The Waterways Experiment Station has made observations 
similar those indicated Fig. 16, both the model and the field.(1) 
However, order avoid too much uplift, 45-degree lip usually 
Efforts use gate with upstream seals and the upstream skin plate 
ing below the bottom members were successful far hydraulic loading 
the gate was concerned. However, the expansion the jet laterally into the 
gate slot intensifies the possibility cavitation the conduit liner. 

For submerged gates with downstream seals, the depth water over the 
gate contributes materially the downpull force. the gate raised and 
the top seal broken, the clearance the downstream side the gate 
significant. The relative clearances the upstream and downstream sides 
can maintain release the column water over the gate, and thus change 


the magnitude the loading. has been found best maintain some water 
the gate well over the gate for more positive control the gate. 

mention was made studies the vibration characteristics the St. 
Lawrence gates. With flat lip, gate tends vibrate vertically 
with the vortex trail shed from the upstream edge. With degree 
the vortices are shed from the downstream edge the gate, thus eliminating 
any areas which pressure pulsations can act. The shape the lips the 
gates Fig. cannot determined, although the basis the tests 
diversion port gates Long Sault Dam, assumed that they are either 
degrees. Other contributing factors the hydraulic loading 
are venting downstream from the gate and entrapment water the 


members. 

Another problem area not mentioned Mr. Bryce’s paper concerns 
emergency closure ports gated areas means stop logs series 
built members. These devices are designed lowered through 
ing water provide barrier flow the event damage gate. 
tests and field experiences indicate that almost impossible install the 
emergency closure devices flowing water. Flow over the top and be- 
neath the stop logs built steel members sets vortex trails which caus 
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the closure devices oscillate violently vertical direction. The ampli- 
tude movement such that the closure device could damaged wedged 
skewed position. 

realized that space requirements prevented the author the basic 
paper from going into detail specific problems. However, hoped that 
his closing remarks the author can expand his tests develop energy 
dissipators and investigate the hydraulic loading gates. 
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GROUND-WATER PROBLEMS NEW YORK AND NEW 
Discussion Biemond 


stated the introduction the paper, there gener- 
ally interchange between surface water and ground water under natural con- 
ditions. There can also interchange artificial processes: surface water 
can brought intermingle with and diffuse ground water artificial 
replenishment, practised Sweden and The Netherlands. The ground 
water captured works this nature mixture natural ground water 
and artificial ground water, which has many respects the same quality 
natural ground water. Artificial recharge accomplished these cases 
spreading shallow ponds canals. The writer does not know examples 
Europe recharge with surface water means wells. 

Intermingling ground water and surface water also brought about 
what called “Uferfiltration” Germany and practised large scale 
along the Ruhr River and the Rhine River. From the paper deduced that, 
the case the Blackstone River near Woonsocket, I., this process, 
which wells for the withdrawal ground water are placed near the river 
bed, being envisaged. 

the purpose the writer report some features resulting from 
this technique, that might interest. 

The alluvial deposits and along river beds—where the river has been ag- 
grading—consist rule considerable depths mixture boulders, 
gravel and coarse sand—a favorable condition for ground-water withdrawal. 
When the wells are placed near the river and reach into the underlying de- 
posits, the water withdrawn would partly percolate out the ground- 
water body adjacent the river bed, but, result the flow 
pattern, mostly through the river bed, and therefore would consist for the most 
part ground water that was recently surface water This 
water would contain the natural impurities river water, such silt and de- 
cayed algae, but also the products artificial pollution, which might or- 
ganic substances, ammonia, iron, and others. Furthermore, there may be, 
result heavy pollution, very low oxygen content. 

Pollution kept away from the Ruhr River much possible; water 
percolating through the gravel and sand the river bed reaches galleries and 
wells satisfactory condition. 

The Rhine River state increasing pollution and the water with- 
drawn from wells bordering the river bed gradually giving more trouble. 

The troubles consist firstly deteriorating quality the extracted water and 
Secondly gradual clogging the river bed. Deterioration quality 
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particularly apparent the taste. Clogging, the opinion the writer, may 
caused high entrance-velocities for the water penetrating into the al- 
luvial deposits, thus introducing insoluble particles into the interstices the 
deposits. might also caused introducing water into 
medium which doubtless contains ferric deposits that might thus into so- 
lution, deposited subsequently interstices where intermingling with 
oxygen-rich water would take place. 

the opinion the writer artificial replenishment, comparison with 
wells located alluvium along rivers, has certain advantages that may 
worthwhile consider. 
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FLOODS THE FLORIDA 


ASCE.—There are few large inland bodies water 
the United States that are subjected hurricane winds that opportunities 
compare their tidal characteristics are rare. Pamlico Sound, eastern 
North Carolina, one the few such waters that can likened Lake 
Okeechobee. This sound land-locked virtue string sand barrier 
reefs, known the Outer Banks, which separate the sound from the Atlantic 
Ocean. has surface area nearly 2,000 square miles, and drains 
area some 25,000 square miles Virginia and North Carolina. Its princi- 
pal outlets the ocean are three relatively narrow and shallow inlets through 
the Outer Banks; thus the normal sound levels can considered inde- 
pendent ocean tides, except for local effects the immediate vicinities 
the inlets. Unlike Lake Okeechobee, Pamlico Sound not “huge saucer”, 
but assumes the form rough parallelogram bounded the north and west 
the relatively steep shores the marshy flatwoods the mainland, and 
the east and south the shoal waters contiguous the Outer Banks. The 
area maximum depth, varying from feet, found narrow reach 
along the northeast-to-southwest major axis the sound. 
Hurricanes are not strangers eastern North Carolina. fact, list 
these storms that have directly affected Pamlico Sound alone since 1900 would 
long the author’s list associated with his Fig. Analysis Fig. 
and its companion list shows that the Okeechobee area crossed two 
general groups hurricane paths. One associated with storms from the 
Atlantic Ocean that generally strike southeastern Florida prior their re- 
curvature. Storms the other group approach from the Caribbean Sea and 
the Gulf Mexico, striking the west side the peninsula after having re- 
curved eastward. noted that the highest wind speeds generally have been 
more closely associated with the storms the former group. Because the 
orientation the Florida peninsula relation its watery environs, ex- 
cessive rainfall not necessarily peculiar either group. However, 
interesting note that the 1924 hurricane produced exceptionally heavy rain- 
fall areas north Lake Okeechobee, whereas the storm center passed 
least miles south the lake; the pronounced eastward swing the path 
this storm (number VI) suggests the proximity stationary front, which 
could account for the widespread area excessive precipitation. 
Because Pamlico Sound situated some eight degrees north and five de- 


grees east Lake Okeechobee, the majority its hurricanes approach from 
the Atlantic Ocean during after recurvature. When hurricane recurves 
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the higher latitudes, its forward speed usually increases considerably, 
carrying the storm rapidly past across the sound; therefore, the regions 
maximum winds not have sufficient time which develop severe wind 
tides over the longer fetches, unless the radius the maximum wind field 
exceptionally large. Three the most severe storms this region have 
amply met this requirement; along the Outer Banks the south side 
Pamlico Sound, the winds and sound tides the “Great Atlantic Hurricane” 
1944 are remembered the most severe half century. 

tropical storm enters its zone recurvature, not unusual for its 
forward speed decrease somewhat; occasion, blocking ridge high 
pressure will cause the storm stall and wander aimlessly. Hurricane Ione 
1955 reached virtual standstill immediately after crossing the North 
Carolina coastline the vicinity Cape Lookout (located about miles due 
south the western end the sound); then proceeded drift slowly and 
erratically westward, covering path distance only miles nine hours 
before veering northeastward and picking speed. During the “stagnant” 
period, the western portion Pamlico Sound experienced winds 100 
miles per hour, and storm tides reached unprecedented heights many muni- 
cipalities that area. 

Not infrequently the varying conditions surface friction presented 
hurricane passing over land masses will cause the eye the storm become 

near the surface, thereby altering surface wind patterns. The scar- 

city surface observations made during the passage Hurricane 

1953 leaves the writer somewhat uncertain the configuration the 
surface eye immediately after crossed the coastline the immediate vicini- 

Cape Lookout. The storm was infant, hurricanes go, its maximum 

steady winds just reaching hurricane force. Weather Bureau summaries 


the storm not suggest the presence unusual characteristics, nor does 
reconstructed path, based available information, reveal significant irregu- 
larities. Had the eye passed undistorted across the western part 
Sound, the maximum setup should have occurred along the north shore. Al- 
though flooding did occur where expected, the highest setup appeared have 
occurred along the northern portion the east shore the sound; the writer 
observed fresh high water line, estimated feet above normal 
level, extending for several miles. Moreover, structural damage, which was 
largely confined area the Outer Banks about miles further north, 

all occurred within the same half-hour period; the damage was inflicted 
southeast wind least miles advance the apparent eye. Curiously 
enough, the path the storm eventually passed miles west the damage 
zone, but the eye was not observed within miles that area. 

Florida also has experienced some unusual behavior the part hurri- 
canes; among the more recent examples were the small ‘but intense Miami 
storm and the looping Cedar Key storm, both The former (number 
XVIII Fig. maintained tight circulation pattern around very small eye 
from Miami Lake Okeechobee; however, the published observations indicate 
considerable reduction wind speed the Okeechobee area, with winds 
only gale force observed the vicinity the lake. After crossing the 
the eye broadened and full hurricane intensity was regained. 

Although the Cedar Key storm September, 1950, was well removed from 
the Everglades, the writer wonders how might have affected the project 
area, assuming the loop, the subsequent path and the associated precipitation 
pattern were superimposed Lake Okeechobee and the Everglades. Cedar 
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Key received 25.2 inches rain little over hours, and depths 
inches were recorded elsewhere along the path the storm. 

The author’s estimates the effectiveness the Central and Southern 
Florida Flood Control Project are impressive, but one wonders whether the 
general public fully appreciates their value. the writer’s contact with 
municipal flooding problems, has observed that having assumed consider- 
able obligation approving bond issue, one derives small comfort from the 
knowledge that the new storm relief system probably will reduce the depth 
water his basement from four feet only two. 
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ROLL WAVES AND SLUG FLOWS INCLINED OPEN 


ESCOFFIER,! ASCE.—The author has made interesting study 
certain waves instability that occur fairly shallow water. Apparently 
either surface tension laminar flow plays essential part the gener- 
ation these waves. the interest providing somewhat more complete 
picture the subject waves instability the writer feels that may 


helpful present his understanding how some waves instability originate. 


The type waves the writer has mind normally occurs relatively 
deep and fully turbulent interesting photograph such roll wave 
having height five feet was published 1936 Holmes(1) together 


with some rather detailed measurements this and other roll waves observed 


him. News received from time time highly destructive roll waves 
that occur canyons the western part the United States. These are 
generally referred “walls water.” 

The writer will use substantially the same reasoning that was used 
1950(2) except that the interest simplicity will assumed that there 
exists uniform distribution velocity the channel cross-section, as- 
sumption that greatly simplifies the results. 

The dynamic equation and the equation continuity, for uniform channel, 
which serve starting point for the derivation are 


and 


where the mean velocity, time, the distance measured down- 
stream, the acceleration gravity, the depth water the lowest 
point the cross-section, the bottom slope, the normal velocity, 
the water-surface width, and the cross-sectional area. 

Eq. (2) now multiplied and the result added Eq. (1) obtain 


The condition for the left-hand side Eq. (3) directional derivative 
the (x,t) plane 
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which yields the following solution for 


(4) 


this back into Eq. leads the equation 


where 
(6) 
and 
(7) 


stage variable which used instead represent the level 
the water the channel. rectangular channel assumes the familiar 
form 


z=viw (8) 
and 
(9) 
are substituted into Eq. (5) there obtained the equation 
(10) 


Use now made the identity 


eliminate the partial derivative from Eq. (10) with the result 
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from which can seen that along the path the (x,t) plane defined 
the expression brackets disappears and 
Fdt 


which the original terminology becomes 


Eqs. (12) are the equations characteristics. The characteristic 
corresponding the upper plus sign represents pulse wave point 
traveling downstream with velocity and shall designated forward 
characteristic. The one corresponding the lower minus sign represents 
pulse wave point traveling with velocity and consequently moving 
either upstream downstream accordingly less than greater than 
shall designated backward characteristic. 

otherwise uniform flow which 


and, therefore 
Along the forward characteristic 


and along the backward characteristic 


Fig. the lines and represent forward characteristics and the 
line backward characteristic. The zone the right assumed 
undisturbed uniform flow. pulse disturbance starts and moves 
downstream along MP. Ordinarily the effect channel friction will re- 
duce the magnitude this pulse moves downstream. However, un- 
usually steep channels the opposite may true and the pulse will continue 
increase magnitude. such case the flow said unstable. 

The line has been taken short enough for the right-hand side Eq. (16) 


Eq. (17) now used eliminate from (15) with the result 
dVo 
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The substitution Eq. (19) into Eq. (18) yields 
= gs dVo 
examination Eq. (20) will show that will increase decrease 
magnitude moves downstream accordingly greater than less 


than one. The requirement for the formation waves instability is, there- 


fore 


(21) 
Eq. (20) can integrated. The result 
(22) 


The same line reasoning with the roles the forward and backward 
characteristics interchanged yields the formula 


(23) 


the requirement for the formation waves instability along the backward 
characteristics. This unusual condition that requires the normal 


Zone disturbed flow 


Zone undisturbed flow 


Fig. 
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velocity diminish the stage increases. However, condition 
that can met closed conduit flowing nearly full. 


remains show that Eq. (21) will upon appropriate substitution yield the 
more familiar formulas general use. rectangular channel infinite 
width 


Therefore 


The formula for the critical slope obtained equating given Eq. (6) 


and 


(26) 


Eliminating and between Eqs. (26) and (27) there obtained 
(28) 


Eq. (21) becomes 


whence 
(29) 
For the Chezy formula and Eq. (29) becomes 
S>4 


which the result obtained Jeffreys(3) and the Manning 
exponent assumed Eq. (29) becomes 


which the result obtained Keulegan and 
also possible derive Vedernikov’s formula from Eq. (21). Veder- 
starts with the equation 
where the hydraulic radius. 


The use this equation would require the substitution the exponent for 
the two exponents the right-hand side Eq. (5). However, this change 
leaves Eq. (21) unaltered. 
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Solving Eq. (30) for there 


Taking the derivative this expression, and then dividing this derivative 
the expression itself yields 


Since 


but 
and therefore 


where the wetted perimeter. 
Letting the velocity pulse moving along forward characteristic 


that 
(33) 


can eliminated from Eq. (32) with the result 


The expression the right-hand side Eq. (34) has been called the Veder- 
nikov number Powell(8) who states that when this number exceeds one, the 
flow ultra-rapid, roll waves form and the flow cannot steady. Veder- 
nikov’s criterion is, therefore, equivalent Eq. (21). However, both have 
been developed here for uniform distribution velocity the channel 
cross-section. Since Vedernikov introduces the velocity-distribution coef- 
ficient into his formula way different from that followed the writer 
his original derivation, the two criteria would differ somewhat for non- 
uniform distribution velocity. 

Craya(9) has shown that Eq. (21) equivalent 


the celerity Seddon, being the normal discharge. Eq. (35) can 
derived substituting 


(35) 
where 
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QodA 


and 


into Eq. (21). the same substitution made into Eq. (23) there obtained 
(37) 


which the corresponding condition for waves instability along the back- 
ward characteristics. 
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REVISED COMPUTATIONS VELOCITY HEAD WEIGHTED 


Discussions Byron Aldridge and Anand 


BYRON ALDRIDGE, ASCE.—The author’s main objection the 
use the present method determining weighted velocity head, 
appears the cumbersome figures obtained. Since only the relative values 
the numerator and denominator are needed, the work can reduced con- 
siderably arbitrarily shifting the decimal point both the area and con- 
veyance columns and making mental calculations The Geologi- 
cal Survey uses this method much its slope-area work. The method 
illustrated below with the example used the authors their presentation. 


Approximate values 


A K 2 rule Kj 
,895 5,000 100 149 1,490 
1,492 


The lines through the and columns are drawn arbitrarily obtain 
small numbers. The decimal can shifted any amount and need not shift- 
the same amount the area column the conveyance column. 
The columns headed “approximate” can obtained rapid mental calcu- 
lation give the location the decimal point for the shifted computations, 
and need not recorded. 
Where there large variation conveyance often more convenient 
use decimals place whole numbers. the above example the decimal 
the area column could shifted places instead which would move 


Proc. Paper 2149, September, 1959, Lara and Schroeder. 


Hydr. Engr., Geological Survey, Surface Water Branch, Tacoma, 
Wash. 
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the decimal one place further the left the column and places further 


the left the remaining columns. 

The value alpha the same either method and correct even though 
the absolute values K3/A2 are not known. 

When the mental calculations show that small areas will not effect 


the last significant figure the total not necessary compute for 


these areas. The data that the writer would record the computation sheet 
for the example shown below. 


Alpha 
1,895 175,000 14.9 
176 3,670 
82.4 935 
14.9 

12.5 


“Velocities the subsections are computed for the final discharge (which 
usually obtained weighting the discharge from subreaches) the 
same manner that shown Lara and Schroeder. 


ANAND,! ASCE.—A study the “revised” formulas recom- 
mended the authors shows that these are essentially the same the “old” 
formula, viz; 


Where the velocity head correction and suffix “p” portion the channel. 
This can proved follows, starting with the following “new” steps the 
author’s. 


Vp P 
(2) col. pp. 


Engr., The Shawinigan Eng. Co., Ltd., Montreal, Quebec. 


4 in 
1 Ca 
4 EAM 


(1) 
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Substituting (4) (3) get 


which the “old” formula given (1). 

for the convenience and saving computational work, the writer fails 
see how working out value from formula (1) one step can more 
cumbersome than filling seventeen column table. velocities and dis- 


charges are required subdivision, they can easily worked out using the 


various values 


Also, objectionable introduce the conveyance term AR2/3 


hydraulic calculations when similar term 1.486 already well 
established. This latter term depends solely channel geometry and has 
predictable shape when plotted for various types channels. Mannings “n” 
can differ from reach reach for the same channel and should used 
separately. 
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DIVISION ACTIVITIES 
HYDRAULICS DIVISION 


Proceedings the American Society Civil Engineers 


NEWS 


December, 1959 


PURPOSE THE HYDRAULICS DIVISION 
(Quoted from the Official Register) 


“The advancement and dissemination knowledge relating the occur- 
rence water nature and its behavior structures, water courses, and 

“In particular the field the Hydraulic Division shall embrace meteorology 
and hydrology they affect the engineer, fluid mechanics engineering 
usage, and applied hydraulics branch engineering science which fur- 
nishes the basis for hydraulic design and for the practical use water the 
different specialized branches hydraulic engineering.” 


EXTRACTS FROM 
MINUTES MEETING THE EXECUTIVE COMMITTEE 
HYDRAULICS DIVISION ASCE 


This meeting the Executive Committee was held Tuesday, June 
1959, the Denver Federal Center, Denver, Colorado. The Committee Chair- 
man, Professor Carl Kindsvater, Those attendance were: 


Professor Carl Kindsvater, Chairman 
Dr. Arthur Ippen, Vice Chairman 

Mr. Harold Martin 

Maurice Dickinson 

Joseph Tiffany, Secretary 


HYDRAULICS DIVISION BUDGET. report received subsequent this 
meeting and dated August 1959 indicates expenditure that date about 
$12,000 the $13,500 authorized for 1958-59. The following budget has been 
submitted the Vice Chairman for 1959-60. 


Note: No. 1959-48 part the copyrighted Journal the Hydraulics Division, 
the American Society Civil Engineers, Vol. 85, 12, December, 
1959. 


Copyright 1959 the American Society Civil Engineers. 
1959-48--1 


1959-48 
| 


1959-48--2 December, 1959 


Executive expenses 3,260 
Technical Committees 16,670 
Total 19,930 


The ASCE Board has not yet acted upon this proposed budget. 


NAME AND PURPOSE THE DIVISION. After considering the various 
names which had been suggested for the present Hydraulics Division Water 
Supply and Hydraulics Division, Hydrology and Hydraulics Division, Hydraulics 
and Water Resources Division, etc. the Committee confirmed its February 
decision retain the present all inclusive name “Hydraulics Division.” 

After carefully considering comments made the Water Resources Co- 
ordinating Committee (WRCC), was decided retain the wording the 
purpose the Hydraulics Division follows: 


“The advancement and the dissemination scientific and engineering 
knowledge all branches hydraulics, hydraulic engineering and water 
resources development. particular this shall embrace meteorology 

and hydrology the sciences dealing with the occurrence water 

the atmosphere, the earth surface and the ground, fluid mechanics 
for the understanding all flow phenomena, applied hydraulics for the 
design and planning hydraulic structures and comprehensive sys- 
tems, and those social, economic and administrative aspects basic 

the conservation and utilization water essential natural resource.” 


The Secretary was instructed prepare letter Society Headquarters 


formally requesting change the stated purpose the Division quoted 


ORGANIZATION WATER RESOURCES COMMITTEES. After careful con- 
sideration the comments the WRCC and the Committee Division 
Activities, was decided proceed immediately with the abolition the 
Committee effective October 1959 and with the establishment 
the four new committee soon they can authorized. The Secretary 


was requested prepare letter Society Headquarters requesting author- 
ization for these actions. 


WATER RESOURCES COORDINATOR, was decided establish position 
Water Resources Coordinator, who would responsible the Executive 
Committee for the activities the four new Water Resources Committees 
until they have been successfully established. would meet with and report 
the Division Executive Committee, and would meet with the new committees 
whenever considers the occasion demands. The Secretary was requested 
prepare letter requesting authority establish this position. 


TECHNICAL COMMITTEE The procedure followed 
technical committees establishing task forces and naming members thereon 
was discussed. Attention technical committee chairmen should called 

the requirement that the establishment new task forces cleared with 
the Executive Committee before any steps are taken name members task 
forces. was further concluded that one individual should asked 
serve more than two technical division committees any one time, and 
that technical committees should check the ASCE Official Register before 
nominating members task forces, insure not selecting individuals already 
several other committees and task forces. 
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Appointments and 1959-1960 membership, administrative committees 
COMMITTEE PUBLICATIONS. Professor Wallace Lansford has agreed 
serve chairman this committee. Mr. James Smallshaw was also se- 
lected for membership this committee. The committee thus will consist of: 


Name Year Termination 
Wallace Lansford, Chairman 1961 
James Smallshaw 1963 
Arthur Ippen 1960 
Harold Martin 1960 


NEWSLETTER Mr. James Ball has agreed serve the 
Newsletter editor for least the ensuing year. 


RESEARCH COMMITTEE. Mr. Jacob Douma was selected the new 
member, replacing Mr. Goodrich whose term has expired. Thus membership 
the committee will be: 


Name Year Termination 
Vito Vanoni, Chairman 1961 
James Robertson 1960 
Clarence Wicker 1962 
Walter Langbein 1963 
Jacob Douma (new member) 1964 


COMMITTEE SESSION PROGRAMS. The committee membership for 
1959-60 will consist the Vice Chairman the Executive Committee, the 
new member the Executive Committee, and the chairman the nine tech- 
nical committees (including the four new water resources committees and ex- 
cluding the Committee Hydrology, abolished), follows: 


Name Year Termination 
Maurice Dickinson 1960 
Eugene Fortson, Jr. 1961 
Arno Lenz 1960 
Emmett Laursen 1960 
Enos Carlson 1960 
Alvin Peterka 1960 
Haywood Dewey, Jr. 1960 


Chairmen the four new Water Resources Committees are selected 
later. (See also this Newsletter, NEW WATER RESOURCES COMMITTEES). 


STANDARDS COMMITTEE, The following terms office were selected for 
the three members this committee previously appointed without definite 
terms office: 


Name Year Termination 
Rolland Carter, Chairman 1961 
Norbert Leupold 1963 
Frank Campbell 1965 
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COMMITTEE THE STEVENS AWARD. Mr. Tate Dalrymple was ap- 
pointed chairman and Professor Victor Streeter was appointed the new 

member replace Professor David Appel, retiring. Membership for the 
next year is: 


Name Year Termination 
Tate Dalrymple, Chairman 1960 
Jacob Douma 1961 
Victor Streeter (new member) 1962 


COMMITTEE THE KARL EMIL HILGARD HYDRAULIC The fol- 
lowing were appointed for two-year concurrent terms this committee: 


Name Year Termination 
William Hiatt, Chairman 1961 
Joseph Howe 1961 
Joseph Bradley 1961 
Appointments and 1959-60 membership, technical committees 


COMMITTEE FLOOD CONTROL 


Name Year Termination 
Arno Lenz, Chairman 1961 
Joseph Perrey (term extended) 1960 
Francis Christian 1962 
Arthur Luecker 1963 
Walter Schulz (new member) 1964 


Joseph Friedkin, resigned. 


COMMITTEE HYDROMECHANICS 


Name Year Termination 
Emmett Laursen, Chairman 1961 
John Ripken 1960 
Norman Brooks 1962 
Frederick Brown 1963 
Douglas Baines (new member) 1964 


Donald Harleman, retiring. 


COMMITTEE HYDRAULIC STRUCTURES 


Name Year Termination 
Alvin Peterka, Chairman 1961 
Glenn Hands 1960 
Harold Pratt 1962 
Rex Elder 1963 
Marvin Webster (new member) 1964 


Jacob Douma, retiring. 
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COMMITTEE SEDIMENTATION 


Name Year Termination 
Enos Carlson, Chairman 1961 
Albert Gildea 1960 
Robert Banks 1962 
Paul Benedict 1963 
Peter Eagleson (new member) 1964 


Alvin Anderson, retiring. 


COMMITTEE TIDAL HYDRAULICS 


Name Year Termination 
Haywood Dewey, Jr., Chairman 1961 
Lindsay Disney 1960 
Irvin Ingerson 1962 
Raymond Boucher 1963 
Charles Bretschneider (new member) 1964 


Eugene Fortson, Jr., retiring. 


Future Plans 
HYDRAULICS DIVISION SESSIONS FUTURE ASCE CONVENTIONS 


New Orleans, La. 7-11 March 1960. Four sessions are planned for this 
Convention. 


Reno, Nevada 20-24 June 1960. Three sessions are planned for this Con- 
vention. 


Boston, Mass. 10-14 October 1960. Six sessions are planned for this Con- 
vention. 


FUTURE HYDRAULICS DIVISION MEETINGS 


Seattle, Washington August 1960. Six sessions are planned for this 
meeting. The exact dates this Conference August 1960 and other ar- 
rangements therefore will handled Dr. Ippen. 

1961. the date this meeting invitation has been received from 
proposed sponsor. 


PRESIDING OFFICERS, FUTURE ASCE CONVENTIONS. The following were 
designated preside Hydraulics Division sessions ASCE Conventions: 


New Orleans Mr. Fortson, Jr. 7-11 March 1960 
Reno Mr. Dickinson 20-24 June 1960 
Boston Dr. Ippen 10-14 October 1960 


Organization Executive Committee 


NEW MEMBER. was noted that Mr. Eugene Fortson, Jr., Chief the 
Hydraulics Division, Waterways Experiment Station, Vicksburg, Mississippi, 


has been appointed for four-year term the Executive Committee effective 
October 1959. 
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ELECTION 1960 OFFICERS. Dr. Arthur Ippen was elected Chairman, 
Mr. Maurice Dickinson was elected Vice Chairman, and Mr. Harold 
Martin was elected Secretary the Executive Committee for 1959-60. 


DUTIES EXECUTIVE COMMITTEE MEMBERS. For the past several 
years, the Vice Chairman has initiated arrangements for technical programs 
ASCE Conventions and Hydraulics Division Conferences for the year ahead 
when will Chairman; during his year Chairman, makes final ar- 
rangements for those same technical programs. Thus, during any given year 
the Chairman has been working programs for the current year and the 
Vice Chairman programs for the year ahead. was decided change this 
arrangement for the year 1959-60, favor the following arrangement. 
The Chairman relieved responsibility for technical programs order 
enable him concentrate the coordination activities the technical 
committees. The Vice Chairman will serve Chairman the Session 
Programs Committee and will make final arrangements for programs for 
the current year. The first-year member the Executive Committee will 
assume responsibility for initial arrangements for programs for 1961. 


EXECUTIVE COMMITTEE CONTACT REPRESENTATIVES COMMITTEES. 


With the beginning the next Society year more flexible arrangement with 
respect contact members will established. This will take into consider- 
ation the availability and proximity Executive Committee members the 
place meeting the technical committees. 


APPRECIATION FOR RETIRING MEMBER. The Secretary was instructed 
record the minutes expression the Division’s appreciation Mr. 
Tiffany for his seven years service the Executive Committee, including 
one year Chairman and three years Secretary. Mr. Tiffany retired 
from the Committee the end this Society year. 


HYDRAULICS DIVISION PLANS FOR RENO CONVENTION 


The Hydraulics Division planning three sessions for the Reno Convention 
held 20-24 June 1960. The Tidal Hydraulics, Hydraulic Structures, and 
Flood Control Committees will each sponsor one session. Messrs. 
Ingerson and Harold Pratt, respectively, are responsible for the programs 
the Tidal Hydraulics and Hydraulic Structures Committees. Information 
not available this time who will organize the session sponsored 
the Flood Control Committee. 


NEW WATER RESOURCES COMMITTEES 


technical committee the Hydraulics Division. This due the fact that 
the Hydraulics Division has assumed responsibilities the field water re- 
sources assigned the Water Resources Report the Horner Com- 
mittee and the Water Resources Coordinating Committee established the 
ASCE Board Direction under the chairmanship Mr. Finley Laverty. 

implementing these responsibilities, the Hydraulics Division Executive 
Committee has asked approval the Committee Division Activities for 
the establishment the following committees: 
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Hydrometeorology 
Surface-Water Hydrology 
Groundwater Hydrology 
Water Resources Planning 


The Water Resources Coordinator, mentioned earlier this Newsletter, 
will work with and coordinate the efforts the new committees during their 
first year operation. 
his HYDRAULICS DIVISION SESSIONS, NEW ORLEANS CONVENTION 
Available information sessions for the New Orleans Convention (7-11 
March 1960) follows: 
Sponsored the Hydromechanics Committee 
VIBRATION HYDRAULIC STRUCTURES, Campbell, Chief, Hydraulic 
Analysis Branch, Army Engineer Waterways Experiment Station, 
Vicksburg, Mississippi. 
THE DAMPING OSCILLATORY WAVES DUE LAMINAR BOUNDARY 
LAYERS, Engleson, Assistant Professor Hydraulic Engineering, 
Massachusetts Institute Technology, Cambridge, Massachusetts. 
FREE SURFACE PHENOMENA TURBULENT FLOW, Dellour, 
Associate Professor Hydraulic Engineering, Purdue University, Lafayette, 
Indiana, and his associate, Mr. McManus. 
Sponsored the Hydrology Committee 
APPRAISAL THE STREAM-GAUGING NETWORK THE UNITED 
STATES, McCall, District Engineer, Geological Survey, Trenton, 
New Jersey. 
tion 
SEEPAGE THROUGH LAYERED ANISOTROPIC POROUS MEDIA, Todd, 
Associate Professor Civil Engineering, and Jacob Bear, Graduate Re- 
search Engineer, University California, Berkeley, California. 


PRESSURE DISTRIBUTION PARTIALLY SATURATED POROUS MEDIA 
DURING STEADY FLOW IMMISCIBLE Scott, Associate 
Professor Irrigation, University California, Davis. 


The Committee Hydraulic Structures will not sponsor session the 
New Orleans Convention; therefore correction page the August 
Newsletter order. 


INTERESTING ASME PROGRAM 


Many ASCE members will interested the Hydraulics Division ASME 
Conference held Houston, 7-9 March 1960. The program will include 
two special symposia, one cavitation highlighting the problems face 
and research needs rather than specific accomplishments. The other 
stall and unsteady flow fluid machinery. 


ASCE 
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FOR YOUR CALENDAR 


ASCE Meetings 


March 7-11, 1960 
June 20-24, 1960 
August 17-19, 1960 
October 10-14, 1960 
April 10-14, 1961 
October 16-20, 1960 
February 1962 

May 1962 

October 15-19, 1962 


February 12-13, 1960 


March 7-9, 1960 
April 4-8, 1960 
April 28-30, 
June 18-21, 1962 


ASCE, New Orleans Convention 

ASCE, Reno Convention 

Hydraulics Division Conference, Seattle 
ASCE, Boston Convention 

ASCE, Phoenix Convention 

ASCE, New York Convention 

ASCE, Houston Convention 

ASCE, Omaha Convention 

ASCE, Detroit Convention 


Non-ASCE Meetings 


Southeastern Regional Meeting, Society 
American Military Engineers, 
Jacksonville, Florida 

Gas-Turbine Power and Hydraulic 
Conference American Society 
Mechanical Engineers, Houston, 
Texas 

Nuclear Congress and Exhibit, 
Coliseum, New York, 

Annual Meeting, American Geophysical 
Union, Washington, 

Fourth National Congress Applied 
Mechanics, University California, 
Berkeley, California 


Deadline dates for Newsletter contributions: January issue November 20; 
March issue Jan. 20, 1960. 


USE THE HYDRAULICS DIVISION NEWSLETTER 


You are urged continue use the Division Newsletter for announce- 
ments, inquiries, personnel news, committee reports, surveys and other items 


NEWSLETTER EDITOR FOR 1960 


Your Newsletter Editor for 1960 Mr. James Ball, Bureau Recla- 
mation, Building 56, Denver Federal Center, Denver 25, The as- 
sistance and cooperation the Hydraulics Division officers, committeemen 
and members furnishing material for the Newsletter 1959 are appreci- 
ated sincerely. 

Best wishes for progressive and HAPPY NEW YEAR. 


GUY ARBUTHNOT, 
Newsletter Editor, 1959 
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